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Abstract 


The A, a and B functions for electron broadening of positive ions, 
that is, with an attractive hyperbolic path, are considered in detail 
for the dipole interaction. Functions for a repulsive hyperbola and 
quadrupole interactions are also discussed. 
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I. Introduction 

When classical path calculations for impact broadening are performed 
a common procedure is to expand the required S-matrix elements to second 
order and to express them in terms of A and B functions. In addi- 
tion, when considering broadening of positive ions by electrons and ions, 
the classical trajectory is a hyperbola. In this report we consider the 
functions corresponding to dipole interactions of electrons in detail. 

Dipole interactions of ions and quadrupole interactions are also mentioned. 

( 3 ) 

For an attractive potential (electron broadening) the functions 
for a dipole interaction are 


A(£,e) + iB(£,e) 


2 r°° rx , 2. , 2 


[ dx[ jy £ e +(e -1) shxshy+chxchy-e (chx+chy) ] ^i£ [ e (shx-shy)- (x-y) ] . (1) 
* ^ (echx-l)^(echy-l)^ 


—00 -00 


For a particle of velocity v and impact parameter p causing a 
dipole transition between levels spaced AE apart of an ion of charge Z 


/l +(p 2 /a 2 ) * a “ Ze2/mv2 


£ ■ ZaAE/iiv 

(notice £ and e are unsymmetrized in equation (2)). 

The real part of equation (1) can be written in terms of modified Bessel 
functions 


( 2 ) 


2 

A(C,e) - (U)V ? [|K> e (?0| 2 +-i£-^-|K u (U)| 2 ] 


( 3 ) 
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The imaginary part cannot be simply expressed and Is calculated 


from the real part by means of the Hilbert transform 


. . 1 f e 

..o - - 

* 


A( 1 4-x| ,e) - A(frf x.e.) 


After the work described in this report was completed, fairly simple 


integral representations for the A and B functions have been found by 
e vi (5) 

Out* voattl fc aikwoo nl ooolir »jrf 1*i4 c ft flirt/ 


S. Klarsfeld. ' Our numerical results agree closely with his B function. 


For completeness, his results are 


#00 

AU,e) - 4C 2 eV 5 dt cos (2£t)ch 2 t K^echt) 

o 


2 2 f 7 *'^ 2EQ 2 

B(C,e) « e d6e^ sin 0 K^esinQ) 

'o 


#00 

- 4C 2 e 2 e 1T? dt sin(2?t)ch 2 t K Q (2?cht) 

J n 


To find the line-width it is also required to integrate the A 


function over impact parameters, to obtain 


*00 

,e) - A«,e') 

• c- 


-e n (Ce)Kj^E)K lc (Ce) 


/ - 2 tt£ 

( - 5e e 


[ dt sin (2£t)sh(2t) K (2£e cht), see reference (5)) . 
in ° 


4 ; 

EX’S 
lsa» 


4 


The corresponding function 


bU 


,e) - | B«, 


-t\ de* 

e > 7 * 


( 6 ) 


( « - i + Ce 2 f d6e 2 ^ 0 sin(20) K (2$e sin $) + 

2 * J 0 0 

#00 

+ U e M dt cos(2€t)sh(2t) K (2£echt), see reference (5)) 
* 


has not been calculated. 


II. Calculation of A(S.e) and a(S,e) 

These are both calculated at the same time in subroutine 
BIGA3(X,PP,A,CAPA) where X ■ £e and PP = £. 

(a) Ce * 10“ 5 , 

A(£,e) - 1 + H 
a(C,e) * - log (Ce) 

(b) 10“ 5 < Ce < 4.0 

this region interpolation in a smoothed table [subroutine 
INTERPA (POUT ,X0UT ,A0UT) ] followed by conversion to A and a [In subroutine 
UNSCRAM (X , P , AOUT) ] is used. 
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The Bessel functions to generete the A's were calculated from a 


program C3 NYU BES4 obtained from SHARE. This routine (called COMBES) 


computes the Bessel function J (x) (with v and x complex) for one 


argument and all orders by using the appropriate recursion relationships 


and normalisation factors. Y^(x) is calculated by summation of J v (x). 


For details see reference 6. 


The formulae relating A and a to the values obtained from the 


smoothed table in INTERPA are 


A(5,e) - «e) 2 e 25 


[*••>' * (*- 'Hi 2 ? 1 * 


a(C,£) - (Ce)e‘^K o (Ce)K 1 (?e)*A0UT(2) . 


(c) £e > 4.0 


In this region e^^Kj^Ce) and e^^K^Ue) are calculated by 


SUBROUTINE KKP(X,P,CAY,CAYP) (x = £e, P = O by use of various asymptotic 


formulae : 


(i) In the range £ 2 ^ 3 (e 2 -l) £ 0.01 


a formula due to Balogh Is used 


.(7) (8) 


K lv (vz) 


exp (• HhL 


,i»q) m y 1 V g) l 

v 4/3 s-0 v 2s ) 


where 5 * v 2 ^ 3 ?, \ c 3 ^ 2 - (z 2 -!) 1 ^ 2 - 


arcsec z. 


K 
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and the coefficients 


*.<«> 


n m Zs-m 

m-0 


C 1 ^ 2 B g (5) 


2s+l , /0 

V V- 

m-0 


a - 1, a 
o 8 


(2s+l) (28+3) • » * (6s-l) 


si (144) 


b - 1, b 

O 3 


6 8+1 


6s-l a s 


and U *1 
o 


lt v r 5 2,., 

V 8 [ 3 v +1] 


.. 2 r 9 . 77 2 385 4, 

U 2* 128 + 192 v * 1152 V 1 


„ 3 r 75 , 4563 2, 17017 4, 85085 6, 

V V [ 1024 + 5120 V + “9216 V + 82944 v ] 


U,“ v 4 [ 4465125 + 94121676 v 2 + 349922430 v 4 +446185740 v 6 + 185910725 v 8 ] 

39813120 


where v - l/(z 2 -l)^ 2 


K lv< vz) " H 


i $ •x-trfe j, ^ 


D g (?) 

..2s" 


where 


s' 1/2c , - 


2s-m+l 


l (-l) E a if 3m/2 V, 

m-0 m 2s “ m 




-7- 


»nd V, 


.3 . 7 2, 

v 1 8 + 24 v J 


V 2 


V 3 


2, 15 . 33 2 , 455 4, 

v 128 + 64 v + 1152 v ] 


(Notice: the negative sign was 

omitted In error in reference 7.) 


v 3 [42525 + 451737 v 2 + 883575 v 4 + 475475 v 6 ; 

414720 


The sums of terms involving the A,B,C and D coefficients are calculated 

in subroutine AB(LSWITCH,NU,ZETA,ZSQ,ACC,BCC,CVC,DVC) . The subroutine is 

originally called with dummy variables so that data can be Initialized (LSWITCH*!) 

Subsequent calls are with LSWITCH-2. Notice in particular that in the sums 

of terms multiplying the Airy functions there is considerable cancellation 

so Double-Precision is used. (There should be at least 16 digits accuracy.) 

The Airy Functions (calculated in AIRY) are obtained directly from tables 

(9) 

given by Abramowitz and Stegun (Table 10.11) which are stored as data, 
and then Interpolated. 

In the region Ce-4.0, comparison of K^(5e) with the COMBES routine 

gave accuracies of better than 5 significant figures, and the accuracy 
improves as (£e) increases. In fact, the Bessel functions were checked 

by direct numerical (Numerov) integration of the defining differential 
equation over an extended range including the region where the COMBES pro- 
gram was used. 

2/3 2 

(ii) In the range C (e -1) < 0.01 a formula due to Erdelyi 


et al. is used 


( 10 ) 


(vz) * j exp^y^ ^(-l)®^ sin t (uH*l)ir/3] r|j«+ 


~ (m+l) 
3 


( 11 ) 


Ii + _L 

2 20 


C 1 = v(z-l) 


6 15 
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C 4 C 2 

fi . fi . JL 

24 24 280 

C 1 C 1 43C 1 
120 + 60 + 4800 


K iv (vz) 


‘ I exp (T)jo ( - 1)0 ainltnH-Dfjrtim + ±) C^fj 

L 


where C* * C C' - 1 

O X 


C£ - v(z-l) 

c 3 - i< c 2 > 2+ n 
c 4 ■ i (c 2 )2+ n c 2 

C 5 ” ^ (C 2 )4+ M (C 2 )2+ 4556' 


-ijstil . rlqj-A)" 

3 -^(f) 3 «> 


Accuracy is in excess of 4 figures even for e - 1 less than 10 . 

(d) If C > 10 - cos ” 1 (l/e)J 

then the following asymptotic form is appropriate 


;,«) - e - C08 ' 1 < 1/e » 


a(£,e) ■ tt/ 2 e 


-2c(/15Ti - cos ” ^(1/e)) 


(13) 


Notice that 5/ e 2 -l » p~ - z and for e large (~.e., p large) the above 
formulae are just the asymptotic formulae for straight line paths 


(i.e., A(z) * 7r|z|e"’ 2 ' z l and a(z) - tr/2 e“ 2 l z ') . 
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III. Calculation of B(£,c) 


The B(C,e) were calculated using the Hilbert transform of equation (4) 
with the A(5,e) functions obtained as in II. The interval (0+,«>) was trans- 
formed into (0+, ASYMVAL) where ASYMVAL is the value where the A(£,e) be- 

—6 

come less than 10 . The Integrand in equation (4) has a discontinuity at 

x * £. The interval (0,£) was thus divided into subintervals (usually 10) 
and integration in each sub-interval used a 10 point Gauss-Legendre routine. 
The interval (5, ASYMVAL) was then also divided into subintervals and 10 
point Gauss-Legendre Integration was used for subintervals of increasing x 
until the integral in the subinterval was less than 10 of the total 
integral. (For example, from £ to 2£» we have steps of 0.25 5 with 10 
point integration, and so on for (say) 4 steps between 2? and 105, etc., 
until the convergence criterion is met.) 

Notice that the integrand does not diverge as x ■+ 0 (for example, it 
tends to -2 as 5 -> 0 in this limit). 

To test the accuracy of the above procedure the double integration for 
the imaginary part of equation (1) was performed numerically with the results 


shown: 

e 

5 

B (Hilbert Transform) 

B (Double Integration) 

1.4 

5.0 

1.849 

1.846 

1.05 

50.0 

3.081 

3.099 

1.01 

10 

-11.84 

-11.89 

1.01 

0.5 

-1.755 

-1.737 

1.01 

0.1 

-0.599 

-0.558 

1.0001 

10.0 

-13.41 

-13.36 

1.0001 

1000.0 

-282.00 

-281.7 

1.0001 

30000.0 

-2170.0 

-2165.5 
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Satisfactory agreement Is seen (better than or two percent), since it is 
not necessary to have the B’s to too great accuracy. 

For a given e, B(£,e) was calculated for increasing K until it 
agreed (to 10 4 ) with the following asymptotic expression^ 2 ^ ^ 

B(£,e) + 5 7 f( 7 r-cos’ 1 (l/e)) + 3| (14) 

2?(e -1) L ^ 2 -l J 

+ - .-£? ■ =- [ (Tr-c° s J_(l /01 ( 9e 6 + i 38e 4 + 152e 2 +16) + (229e 4 +592e 2 +124) 1 

25 3 (e 2 -l ) 5 [ 8^2-1 24 J 

2 . 1/2 

Notice that for e large this tends to Tr/4£(e -1) (* tt/4z) which is 

the straight line result. 

After B(£,e) was calculated and tabulated, the data were smoothed. 
The smoothed data are used in subroutine BFINAL. The smoothed values which 
are stored in the tables will here be referred to as and B,j. Then: 

(a) (i) 1.0001 £ e < 1.3 


B(£,e) = 


[- 


ft,-g>«- 1 a/0) 12 — ■ - i + 3 
/cEi . 


[0r/2)(e 2 -l) 1/2 ] 


B T (5 T ,t) 


(15) 


If > 12.5 the asymptotic value of equation (14) is used. For 12.5 

interpolation in the table is used with linear interpolation in log (l/(e-l)) 
for e and quadratic interpolation for 
(ii) e < 1.0001 

This is an unimportant region, and use of the value of B T (£,j,,e) at 
e ■ 1.0001 is sufficient. Notice however, that when e = 1, B diverges as 
as -r 2 (2/3) ~ 2 j 2 C 2 ^ 3 * - 2.88 £j 2 ^ 3 ( see re f erence 5). 


* ■«'***M f 
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(b) Uses 2.0 

No smoothing is used in this region and B(£,e) is stored directly. Quad- 
ratic interpolation in the 5 and e directions are used. For K > 30 
the asymptotic value of equation (14) is employed. 

(c) 2.0 < e < 00 


(ir-cos” 1 (l/e))-^~=)' + 3 

r»(e 2 -l) 3 l 


B (5,e) * 


f -1/1/ ns( 2+£2 ) f ^ 

(tt-cos (l/e))*^==r" + J 

L /e*-l _ 

r,u 2 -D 3/2 ' 

4e 2 


b 2 (c/e -i)b t (c t ,o 


(16) 


The function B 2 (C^l) is the straight line B-f unction. 

If (i) £ t < 0.0025, we put B(p,e) = 0 

(ii) 0.0025 < 5 £ 5.0 and e < 100, table interpolation in 

takes place with quadratic interpolation in the f ; ^ direction and 
linear interpolation in log e for e. For e > 100 in this range of 


B(C,e) = B 2 (c/e 2 -l) 


(ill) > 5.0. The asymptotic value for B(£,e) of equation (14) 


is used. 


In the important regions for the function B(C,e), K and e are usually 


such that B(£,e) is positive. The regions of negative B(£,e) (see, for 
example, the figure in reference 11) is not very important. Also, when e 


► 1 * 
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gets small the asymptotic formula (equation (14)) is usually sufficient and 
for £ small, e is often sufficiently large that the straight line results 
are valid. 

In the region of positive B(£,e) the accuracy should exceed 2%; in 
the negative region the accuracy is expected to decrease as 5 -*■ 0 and 
e -► 1 (however, comparison with reference 5 shows good agreement even 
when e - 1.001). 


IV . Dipole Functions for a Repulsive Interaction 
The functions A^^e) and a R (C,e) are very easily obtained from the 


attractive case, 


(3) 


since 


A R (?,e) = e~ 27r? A(£,e) 

and 

a R (C,e) * e -27T? a(£,£) . 

The B function for the repulsive case has not been calculated, since inelastic 

collisions with ions are negligible and £ is always very large. The 

(3) 

asymptotic form is therefore sufficient, thus 


B R a,e) 


2£(e 2 -1) 2 



-1 

COS 



(notice: there is a typographical error in reference 3). 

This result does not diverge as e 1. 


V. Quadrupole Broadening Functions 

This section is included for the sake of completeness. It is found ^ 

that quadrupole terms are only important for levels which are essentially 

degenerate, thus we need only to evaluate the function corresponding to 

A^(z)^ in the limit of £ = 0 (i.e. , AE = 0) . 

(3) 

Thus, for an attractive potential 
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(S.e) 1(1 f (7r-cos“ 1 (l/e)) l 2 ) 

P - f 


and 


V t,e) (sf * " h(e) 


■Mi 


(ir-coe^d/e)) 2 

/ 2 , x 2 
(e -1) 


- (tt— cos 


1 (l/e)) ^ j -cos” 1 (l/e) y- 


i \ 
l72 ‘ (c 2 -1) 3/2 y 


+ 1 + i ( x -cos (1/e)) 2 + 

2(e -1) * Z 3e 

Notice that In the limit of e large, both functions reduce to the straight 

line results (i.e., (0) * 4/3 and a^ ■ y (^~^) 2 ) 

For the repulsive potential 

Hs<e) = 7 |^7 t ' (Xd <2 / 2> ] + I7 } 


and 


/ % 1 [l (cos”* (1/e)) . “I/-. / \/ tt -1/, /.x 1 

(O— yhr V i +cos ( 1 / e )( y " cos w e ) " 


*R W/ 2a 2 L 2 


/ 2 « v2 

(e -1) 


<e 2 -l) 1/2 


d? 71 ) 


+ j ( j - cos”*(l/e)) 2 + 3c 


? ] • 


2 (e -1) 

Again these tend to the expected limit for e large. The limit as e 
(i.e., as p -► 0) is more interesting. In particular, 


-* 1 


h R (e-*l) 


Hi - *] 


which is consistent with the ^2^ function of Alder et al. (reference 4 


pg. II. E 72). 
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Yet, 


h(e-KL) 


2 2 

it a 

/ 2 

4p 


which diverges as p -*■ 0. This is a point worth noting, since Cooper and 
Oertel^ 12 ^ used the relationship (2.38) of Biedenharn and Brussard, 
which implies 


f /M m I ^ I f /t\ 

r E2 ^'Attractive £2^ ^'Repulsive 

Since the f_, 0 functions are equivalent to putting e * 1, it is obvious that 

Biedenharn and Brussard's relation is in error for quadrupole interactions, 

although it holds for dipole interactions (i.e., f £ ^ Attrac ti ve * 


El Repulsive^* 

The reason that f E1 Attractive does not diverge as p -> 0 is that 

the hyperbola is strongly curved, and the axial component of the electric 

field changes in direction during the collision and cancellation occurs. 

1 2 

Whereas the quadrupole interaction (proportional to ^Ocos 0-1)) does not 
change sign and there is no cancellation to prevent divergence in the 
attractive case (as 0 goes tt 0 -*■ —it) . There is no divergence in the 
repulsive case since the distance of closest approach is 2a. 
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FORTRAN PROGRAMS FOR a(5,e), A(C,e), B(C,e) 

All page references In the following sub-programs are with respect to 
this report. If the input parameters to these routines are outside the 
Intended range for this report, error messages will be printed, (for example, 
e should be greater than unity) . 

Programs for a(£,e), A(£,e) are listed on pages 18 ff. (see page 17). 
Programs for B(£,e) are listed on pages 34 ff. (see page 33). 


a(€,e), A(£,e) are calculated by a FORTRAN sub program SUBROUTINE BIGA3. 

For example: 

CALL BIGA3 (X , PP , A , CAPA) 
where: x ■ £e 

PP - £ 

A - a(£,e) 

CAPA « A (£,e) 

SUBROUTINE BIGA3(X,PP,A,CAPA) requires the following routines: 

SUBROUTINE INTERPA(POUT,XOUT,AOUT) 19 

SUBROUTINE UNSCRAM (X,P,AT) 23 

REAL FUNCTION KO(X) & REAL FUNCTION Kl(X) 24 

SUBROUTINE KKP(X,P,CAY,CAYP) 25 

SUBROUTINE KIDKI(NU,Z,CAY,CAYP) 26 

SUBROUTINE AB(LSWITCH,NU,ZETA,ZSQ,ACC,BCC,CVC,DVC) 28 

SUBROUTINE AIRY(X,ERI,ERIP) 30 

FUNCTION FACT (N) 32 


non non non non non 
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SUBROUTINE BIGA3(X*PP.A*CAPA) 001 

002 

10 - 25 - 1970 003 

004 

DIMENSION AOUT ( 2 ) 005 

COMMON / SUB/ PI*HALFPI 006 

DATA( PI * 3*141592654) *(UP*l0*0) 007 

P « ABS (PP) 008 

EPSILON ■ X/P 009 

A ■ 0*0 010 

CAPA * 0*0 Oil 

IF ( P *GT • X) RETURN 012 

IF ( ABS (EPSILON-1. 0) .LT.1.0E-50) GO TO 100 013 

ACOSINE * ACOS ( 1.0/EPSILON) 014 

ROOT = SORT (EPSILON**2-1.0> 015 

ASYMVAL * UP / (ROOT - ACOSINE) 016 

IF ( X *GT. 1.0E-5) GO TO 50 017 

018 

P * EPSILON LE 1.0E-5 SEE PAGE 4, II (A) 019 

020 

40 CAPA « 1*0 + PI * P 021 

A * - ALOG(X) 022 

RETURN 023 

50 IF ( P *LT. ASYMVAL) GO TO 100 024 

025 

SEE PAGE 8. II (D) 026 

027 

TEXP a EXP (-2.0*p * (ROOT - ACOSINE)) 028 

CAPA s PI * P * ROOT * TFXP 029 

A = (PI/2.0> * TEXP 030 

RETURN 031 

100 IF(X.GT.4.0> GO TO 125 032 

033 

1.0E-5 LT P * EPSILON LE 4 SEE PAGE 4* 1 1 (B) 034 

035 

CALL INTERPA(P»X*A0UT) 036 

CALL UNSCRAM(X*P»AOUT) 037 

CAPA = AOUT ( 1 ) 038 

A = AOUT ( 2 ) 039 

RETURN 040 

041 

P * EPSILON GT 4 SEE PAGE 5, 1 1(C) 042 

043 

125 CALL KKP(X*P*CAY*CAYP) 044 

XSQ a x**2 045 

COEF =1*0 046 

CAPA = XSQ * COEF *(CAYP**2 + ((XSQ - P**2)/XSQ) * CAY**2) 047 

A = (-COEF) * X * CAY * CAYP 048 

RETURN 049 

END 050 
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SUBROUTINE INTERPA(POUT* 

XOUT.AOUT) 




051 

DIMENSION A( 29*2 

9,2 ) , AOUT ( 2 ) »T ( 2 

,2) ,P(29) ,X(29> 



052 

EQUIVALENCE! P*X) 







053 

DATA! (A! I ) , 1 = 1,377) = 






054 

# 

1*02430 

* 28(0.0)* 






055 

ft 

1*01591 

,1.04513 

* 27 { 0*0 ) t 





056 

# 

1*01609 

♦1.03190 

,1.05892 

, 26(0.0) , 




057 

# 

1*01686 

,1.02995 

,1.04483 

,1.06904 

, 25(0. 

0 ) , 



058 

* 

1.01759 

*1.03041 

,1.04160 

,1.05544 

,1.07693 

, 24(0. 0>* 


059 

# 

1*01804 

,1.03129 

,1,04152 

,1.05150 

,1.06421 

,1.08321 

, 23(0.0), 

060 

ft 

1.01844 

,1.03232 

,1.04262 

,1.05129 

, 1.06029 

,1.07183 

,1.08862 

• 22(0.0), 

061 

ft 

1.01879 

,1.03334 

*1.04412 

,1.05257 

, 1.06010 

,1.06823 

,1.07865 

,1.09351, 

062 

ft 

21 <0.0> 

, 







063 

ft 

1.01908 

,1.03426 

,1.04569 

,1.05444 

,1.06159 

,1.06824 

,1.07557 

*1.08495, 

064 

ft 

1*09810 

, 20(0 

.0) , 






065 

ft 

1.01933 

,1.03510 

,1.04722 

,1.05651 

,1.06382 

,1.06999 

,1.07589 

,1.08249, 

066 

ft 

1.09091 

,1.10258 

, 19(0.0 

) , 





067 

ft 

1.02066 

,1.03996 

,1.05707 

,1.07213 

,1.08528 

,1.09665 

,1.10640 

,1.11467, 

068 

ft 

1.12161 

,1.12736 

,1.15700 

, 18(0. 

0) , 




069 

ft 

1.02112 

,1.04189 

,1.06132 

,1.07945 

» 1.09631 

♦1.11192 

,1.12632 

,1.13957, 

070 

ft 

1.15171 

,1.16277 

,1.22534 

,1.24052 

, 17(0.0 

) , 



071 

ft 

1.02131 

,] .04285 

, 1 • 063 5 1 

,1.08333 

, 1.10227 

,1.12038 

,1.13764 

,1.15407, 

072 

ft 

1.16966 

,1.18444 

,1.29128 

,1.33813 

,1.34792 

, 16(0. 

on 


073 

DATA! (A 

( I ) ,1=378,580) = 






074 

ft 

1.02139 

,1.04338 

,1.06480 

,1.08564 

,1. 10591 

,1.12558 

♦1.14467 

,1.16317, 

075 

ft 

1.18107 

,1.19836 

,1.33921 

,1.42630 

,1.46751 

,1.46986 

, 15(0. 

0) , 

076 

ft 

1.02140 

,1.04368 

,1.06561 

,1.08713 

, 1.10829 

,1.12904 

,1.14937 

,1.16929, 

077 

ft 

1.18881 

,1.20788 

,1.37430 

,1.49617 

,1.57423 

,1.60869 

*1.59891 

,14(0.0) , 

078 

ft 

1.02140 

,1.04388 

, 1.06615 

,1.08816 

, 1.10995 

,1.13146 

,1.15270 

,1.17367, 

079 

ft 

1.19434 

,1.21473 

,1.40075 

,1.55146 

,1.66349 

,1.73305 

,1.75649 

,1.73012, 

080 

ft 

13(0.0 > 

, 







081 

ft 

1.02136 

,1.04400 

,1.06652 

,1.08891 

,1.11112 

,1.13323 

,1.15516 

,1.17691, 

082 

ft 

1.19847 

,1.21983 

,1.42125 

,1,59590 

♦1.73797 

,1.84110 

,1.89899 

,1.90689, 

083 

ft 

1.86049 

, 12(0*0 

) , 






084 

ft 

1.02131 

, 1.04406 

,1.06679 

,1.08945 

,1.11205 

,1.13455 

,1.15702 

,1.17937, 

085 

ft 

1.20164 

,1.22377 

,1.43756 

,1.63222 

, 1.80065 

,1.93481 

,2.02644 

,2 *06868 , 

086 

ft 

2.05718 

,1.98837 

, 11(0.0 

) , 





087 

ft 

1.02120 

,1.04410 

,1.06692 

,1.08986 

,1.11269 

,1.13558 

,1*15848 

,1.18132, 

088 

ft 

1.20416 

,1.22690 

,1.45079 

,1.66240 

,1.85394 

♦2.01643 

,2.14026 

,2.21674, 

089 

ft 

2.23961 

,2.20571 

,2.11291 

, 10(0.0 

) , 




090 

ft 

1.02117 

,1.04406 

,1.06707 

,1.09012 

,1.11322 

,1.13645 

,1.15966 

,1.18285, 

091 

ft 

1.20615 

,1.22942 

,1.46173 

,1.68781 

,1.89969 

,2.08797 

,2.24211 

,2.35209, 

092 

ft 

2.40974 

,2.41009 

,2.35152 

,2.23381 

, 9(0.0) 

) 



093 

DATA! (A 

( I ) , 1=581 ,725) = 






094 

ft 

1.02105 

,1.04411 

,1.06716 

,1.09042 

,1.11371 

,1.13709 

,1.16052 

,1.18414, 

095 

ft 

1.20777 

,1.23146 

,1.47089 

,1.70949 

,1.93937 

,2.15104 

,2.33355 

,2.47584, 

096 

ft 

2.56824 

,2.60378 

*2.57901 

,2.49413 

,2.35102 

, 8(0.0 

) • 


097 

ft 

1.02094 

,1.04406 

,1.06722 

,1.09059 

,1.11399 

,1.13758 

•1.16132 

,1.18519, 

098 

ft 

1.20915 

,1.23329 

,1.47866 

,1.72815 

,1.97409 

,2.20704 

,2.41595 

,2.58915, 

099 

ft 

2.71579 

,2.78714 

,2.79768 

,2.74572 

,2.63333 

,2.46466 

, 7(0.0 

) , 

100 

ft 

1.02091 

,1.04401 

,1.06724 

,1.09066 

,1.11431 

,1.13802 

,1.16203 

•1.18615, 

101 

ft 

1.21021 

,1.23460 

,1.48537 

,1.74442 

,2.00469 

,2.25706 

,2.49051 

,2.69309, 

102 

ft 

2.85310 

,2.96041 

,3.00755 

,2.99064 

,2.90982 

,2.76910 

,2.57492 

, 6(0.0), 

103 

ft 

1.02Q83 

,1.04411 

,1.06715 

,1.09083 

,1.11459 

,1.13856 

,1.16256 

,1.18669, 

104 


■V Wav 
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* 

1.21134. 

1.23601 

,1.49111 

,1.75071 

,2.03181 

,2.30191, 2. 55825, 2. 78865, 

105 

# 

2.98096. 

3.12392 

*3.20844 

,3.22855,3.18212,3.07112,2.90149,2.68203* 

106 

# 

5(0*0) • 






107 

# 

1.02090. 

1,04378 

.1.06724 

,1.09104 

,1.11454 

, 1 .13890 *1.16303,1.18758, 

108 

* 

1.21196, 

1,23713 

,1.49615 

,1.77124 

,2.05602 

,2.34241,2.61993,2.87672, 

109 

# 

3.10014. 

3.27818 

,3.40037 

.3.45893 

,3.44950 

,3.37175,3.22954,3.03065, 

110 

# 

2.78622, 

4(0.0) 

) 




111 

DAT A ( <A< 

I ) , 1*726,841)* 




112 

# 

1.02099. 

1.04371 

,1.06725 

,1.09068 

,1.11514 

,1.13909,1.16363,1.18793, 

113 

* 

1.21296, 

1.23792 

,1.50071 

,1.76256 

,2.07778 

,2.37914,2.67645,2.95810, 

114 

# 

3.21137, 

3.42368 

,3.58350 

,3.68142 

,3.71109 

,3.66989,3.55933,3.38507, 

115 

# 

3.15673, 

2.88771 

, 3(0.0) 

, 



116 

# 

1 • 02o57 ♦ 

1.04382 

,1.06761 

,1.09076 

,1.11500 

,1.13916*1.16393,1.18825, 

117 

# 

1.21354, 

1.23878 

,1.50484 

,1.79252 

,2.09750,2.41269,2.72836,3.03342, 

118 

# 

3.31535, 

3.56097 

,3.75807 

,3.89584 

,3.96617 

,3.96430*3.88937,3.74471, 

119 

# 

3.53777, 

3.27990 

,2.98672 

» 2(0.0), 


120 

# 

1.02129, 

1.04348 

,1.06793 

*1.09154 

,1.11553 

,1.13967,1.16498,1.18886, 

121 

# 

1.21459, 

1.23967 

,1.50910 

,1.80229 

,2.11528 

♦2.44294,2.77604,3.10333, 

122 

# 

3.41255, 

3.69055 

,3.92435 

,4. 10219 

,4.21429 

,4.25396,4.21808,4.10766, 

123 

# 

3.92782, 

3.68770 

,3.40033 

,3.08343 

, 0.0, 


124 

# 

1.01947, 

1.04329 

,1.06764 

♦1.09038 

,1.11519 

,1.13990 *1.16435, 1.1 8840, 

125 

# 

1.21504, 

1.23930 

,1.51168 

,1 .80999,2.13168,2.47112*2.82010,3.16834, 

126 

# 

3.50368, 

3.81294 

,4.08281 

♦4.30056 

,4.45513 

,4.53806,4.54414,4.47206, 

127 

* 

4.32459, 

4.10864 

♦3.83497 

♦3.51818 

,3.17802) 

128 

DATA ( (A( 

00 

!) 

2,1218)= 




129 

* 

1.03 551, 

20(0.0 

) , 




130 

* 

1.01816, 

1.03052 

, 27(0. 

0) * 



131 

# 

1.01928, 

1.03377 

,1.04236 

, 26(0.0), 


132 

# 

1.01993, 

1.03578 

, 1 * 04666 

1.05235 

, 25(0.0), 

133 

# 

1.02036, 

1.03718 

,1.04968 

1.05769 

1.06103 

24(0.0) , 

134 

* 

1.02060, 

1.03813 

,1.05187 

1.06163 

1.06725 

1.06860, 23(0.0), 

135 

# 

1.02077, 

1.03888 

, 1 . 05360 

1.06474 

1.07217 

1.07574,1.07538* 22(0.0), 

136 

* 

1.02092, 

1.03949 

,1.05500 

1.06727 

1.07617 

1.08158,1.08338*1.08150, 

137 

# 

21 (0.0) , 






138 

# 

1.02104, 

1.03999 

,1.05616 

1.06938 

1.07952 

1.08644,1.09006,1.09028, 

139 

# 

1.08706, 

20(0.0 

) , 




140 

# 

1.02114, 

1.04042 

,1.05715 

1.07H7 

1.08235 

1.09058,1.09575,1.09776, 

141 

# 

1.09658, 

1.09213 

, 19(0.0 




142 

* 

1.02164, 

1.04268 

,1.06241 

1.08074 

1.09758 

1.11284,1.12645,1.13831, 

143 

# 

1.14837, 

1.15655 

,1.12625 

18(0. 

) * 


144 

# 

1.02181, 

1.04357 

,1.06456 

1.08470 

1.10392 

1.12216,1.13935,1.15545, 

145 

# 

1.17039, 

1.18410 

,1.24359 

• 14468 

17(0.0 

, 

146 

* 

1.02186, 

1.04401 

,1.06570 

1.08686 

1 • 1 0*740 

1.12731,1.14654,1.16502, 

147 


1.18270, 

1.19955 

,1.31223 

1.30094 

1.15601 

16(0.0) ) 

148 

DATA ( ( A ( 

I ) , 1=1219,1421 ) 




149 

* 

1.02184, 

1.04428 

, 1.06640 

1.08821 

1.10960 

1.13057,1.15110,1.17110, 

150 

# 

1.19056, 

1.20945 

,1.35741 

1.40816 

1.34222 

1.16351, 15(0.0), 

151 

# 

1.02184, 

1.04440 

,1.06689 

1.08908 

1.11109 

1.13283,1.15425,1.17531, 

152 

# 

1.19601, 

1.21630 

,1.38941 

1.48627 

1.48351 

1.37362,1.16872, 14(0.0), 

153 

# 

1.02178, 

1.04447 

,1.06713 

1.08976 

1.11219 

1.13442,1.15651,1.17837, 

154 

* 

1.20000, 

1.22136 

,1.41323 

1.54568 

1.59416 

1.54474, 1.39845,1.17245, 

155 

* 

13(0.0) , 






156 

# 

1.02166, 

1.04449 

,1.06739 

1.09018 

1.11294 

1.13561,1.15824,1.18076, 

157 

* 

1.20301, 

1.22522 

,1.43171 

1.59239 

1.68307 

1.68633,1.59580,1.41866, 

158 
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* 1*17519* 12(0.0) * 159 

* 1*02166* 1.04453 *1.06754 *1.09052 » 1 • 1 1364 » 1* 13655 * 1 . 15957 *1.18253* 160 

* 1*20539, 1.22827* 1.44641 *1.63001* 1.75601 *1*00515, 1.76629 *1.63919* 161 

« 1 *43546*1. 17724 » 11(0.0)* 16? 

* 1.02141,1.04453,1.06760,1.09073,1.11401*1.13724,1.16073,1.10405, 163 

* 1.20726,1.23061 *1.45829, 1.66103, 1.81689 *1.90610 *1.91463 *1.83659* 164 

* 1.67668,1.44971,1.17878, 10(0.0)* 165 

* 1.02154,1.04466,1.06773,1.09101*1.11431,1.13781*1,16167,1.18500* 166 

* 1.20900*1.23277,1.46821 ,1.6869 1,1 >,86842 ,1.99298. 2. Q4455, 2. 01356, 167 

* 1.89904,1.70948,1.46196,1.17995, 9(0.0)) 168 

DATA( ( A( I ) * 1*1422 * 1566 ) ■“ 169 

* 1.02143,1.04419,1.06799,1.09134,1.11470*1.13845,1.16238,1.18628, 170 

* 1.20998,1.23427,1.47670*1.70904,1.91258,2.06838,2.15919,2.17266, 171 

* 2. 10360,1.95503,1.73848,1.47263,1.18084* 8(0.0), 172 

* 1.02177,1.04453,1.06771 ,1.09109,1.11448,1.13855*1.16303,1.18692, 173 

* 1.2 116 1,1. 2 3539, 1.483 7 1,1. 72775, 1.95107, 2. 13454 *2. 26100 ,2* 3 1620* 174 

* 2.29174, 2.18602, 2, 00 561,1.76437, 1,4820 1*1.181 52* 7(0.0), 175 

* 1.02187, 1,0436 2 ,l,0677o* 1.09111,1, 11 49 1*1. 13893*1. 16296 * 1 .18799* 176 

* 1.21261,1.23668*1.48957*1.74436,1.98442*2.19275*2.35181*2.44624, 177 

* 2.46485,2.40272*2.26190*2.05163,1.78767,1.49034*1.18204*6(0.0) * 178 

* 1.02162,1.04332,1.06811,1.09138,1.11511,1,13905,1.16301*1.18855, 179 

* 1.2 1363, 1,238 12, 1,49560 *1.75893, 2. 01412, 2. 24457 ,2. 43343, 2. 56431* 180 

* 2.62444,2.60575,2.50661,2.33203,2.09374,1.80878*1.49779*1.18241* 181 

* 5(0.0) » 182 

* 1.02124, 1.04397,1.06652 ,1,09179 ,1.1 1637, 1.14017, 1.16308, 1.18764, 183 

* 1,21348,1.23787,1.50034,1.77104,2,04050,2.29095,2.50707,2.67206, 184 

* 2. 77175, 2. 79594, 2 ,73962 *2. 60412, 2. 39718, 2. 13246, 1,82803 ,1.50451, 185 

* 1.18269, 4(0.0)) 186 

DATA ( (A( I ) ,1=1567,1682)® 187 

* 1.01929 ,1.04 363 ,1.06923 ,1.09123, 1.1 1424, 1.13799, 1.16626*1. 19056, 188 

* 1.2 149 1,1. 23910, 1.50470, 1,78326, 2.0640 1,2.33275,2 .57374, 2 .77047, 189 

* 2. 90807, 2. 97400, 2. 96 105, 2. 86693. 2 ,69585 ,2. 45789, 2. 16823, 1.84567, 190 

* 1.51059,1.18288,3(0.0), 191 

* 1.02469,1.04373,1.06682,1.09332,1.11610,1.14173,1.16362,1.18787, 192 

* 1.21402,1.24167,1.50610,1.79211,2.08525.2.37016,2.63485,2,86092, 193 

* 3.03431,3.14110,3.17124,3.12000,2.98824,2.78240,2.51467,2.20144, 194 

* 1.86191,1.51614,1.18301, 2(0,0), 195 

I * 1.02591,1.04106,1.06408,1.09392,1.11979,1.14194,1.16971,1.19354, 196 

* 1.21370, 1.23849, 1.50943, 1.80317, 2. 10493, 2. 40467,2. 69049,2.94453, 197 

* 3. 15189, 3. 29750 *3. 37061, 3. 36323, 3. 27330, 3. 10400, 2. 86430, 2. 56794, 198 

* 2.23234,1.87694,1.52122,1.18309, 0.0* 199 

* 1.01726,1.04364*1.06638*1.08573*1.11662*1.14342*1.16640*1.18582* 200 

* 1.21 444, 1,23 900*1.51591,1.81069, 2. 11952*2.43593*2 .74149 *3.02169* 201 

* 3.26099*3.44460*3.55961 *3.59674*3.55037*3.42113*3.21465*2.94187* 2 02 

* 2.6180A *2.26124*1. 8 90 88*1.52590 *1.18312) 203 

DAT A(P S .02*.04,.06».08».1» .12, ,14*. 16, .18, ,20, .4*. 6*. 8*1. 0*1. 20* 204 

1 1.40*1.60*1. 8o*2. 0*2. 20,2.40*2.60*2. 80*3.0*3.20*3 .40 *3.60* 205 

2 3.80*4.00) 206 

M = N = 29 207 

DO 10 I * 1 *M 208 

IF ( POUT .LE.P ( I ) ) GO TO 15 2J9 

10 CONTINUE 210 

GO TO 3000 211 

15 IP « I 212 


~ * 


n n n non nnn nnnnnn 
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DO 20 I * 1 *N 
IF(XOUT.LE.X< I>) 60 TO 25 
20 CONTINUE 
GO TO 3000 
25 IX « I 

IF ( IX. EQ. IP) GO TO 200 
H ■ X ( I X ) - X ( I X— 1 ) 

PP « (XOUT-XUX-in / H 

IX .NE» IP 

IF { IP.GT.l) GO TO 100 

IP » 1 0 .LE. P .LE. .02 X(IX-l) .LE. X .LE. X(IX) 4 PTS 

C = P C I P ) 

QO » POUT / C 
TEMPP * 1.0 " PP 
TEMPO * 1.0 - QQ 
TA a TEMPP * TEMPO 
TB * PP * TEMPO 
TC * QQ * TEMPP 
TD * PP * QQ 

F00 * 1.0 F10 = 1.0 

DO 30 I * 1*2 

AOUT(I) = TA + TB ♦ T0*A (I P . IX-1 » I ) + TD*A< IP# IX* I ) 

30 CONTINUE 
RETURN 

REMAINDER OF TABLE - 4 PTS. - NO PT. ON DIAGONAL 

100 C « P( IP) - P( IP-1 ) 

QQ = (POUT-PI IP-1) ) / C 
TEMPP * 1.0 - PP 
TEMPO » 1.0 - QQ 
TA a TEMPP * TEMPO 
TB « pp # TEMPO 
TC a QQ » TEMPP 
TD a PP # QQ 
DO 130 I = 1*2 

AOUT ( I ) a TA * A( iP-l.IX-lf I ) ♦ TB * A<lP-l*IX.I ♦ TC # A( IP* IX-1 • I 
1) + TD # A( I P * I X * I ) 

130 CONTINUE 
RETURN 

IX = IP = 1 o .LE. X AND p .LE. 0.02 ? PT 

200 IF ( I X .NE. 1) GO TO 250 
QQ * POUT/ .02 
PP * XOUT/.02 
TA * 1.0 “ PP - OQ 
DO 210 I = 1*2 


213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 
22 5 
226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 


£36 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 
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AOUT(I) • TA + PP ♦ QQ * A(IXtIP*I) 267 

210 CONTINUE 268 

RETURN 269 

270 

IX * IP 271 

272 

250 H = X(IX) - XIIX-11 273 

C * PC IP) - PC IP-1) 274 

PP » C XOUT-X C IX-1 ) ) /H 275 

00 » (POUT-PC IP-1) ) / C 276 

TA • 1.0 - PP - 00 277 

00 2701 * 1*2 278 

AOUTCI) * TA * AC IP— 1* IX*! ) + PP * A C IP-1 * IX-1 * I ) + OQ * A(lP*IX*l 2 79 
1 ) 280 
270 CONTINUE 281 

RETURN 282 

3000 PRINT 650 *POUT *XOUT ♦ IP » I X 283 

650 FORMAT C * P*X»* 2F20.5t2I10> 284 

CALL EXIT 285 

END 286 


SUBROUTINE UNSCRAM CX*P*AT> 287 

288 

CONVERSION TO A* CAP A 289 

290 

REAL KO*Kl 291 

DIMENSION AT C 2 ) 292 

XSO = X**2 293 

FKO = KOI X) 294 

C K1 * K1 ( X ) 295 

if,"* 5 = EXP(2.0*P> 296 

EPS"* = XSQ/P**2 297 

UA » l.o-P**2/C CP**2+1.0)*EPSO) 298 

EPSO; = EPSQ 299 

UB * :EPSQ2 - 1.0 + (1.0/CP +1.0) ) ) / EPSQ2 300 

TCA * XSQ * TEMP * (FK1**2 + UA * FKo**2) * UB 301 

,'A * X * TEMP * FKO * FK1 ?02 

AT ( 1 ! « AT ( 1 ) * TCA 303 

AT C 7 ) = ATC 2 ) * TA 304 

RET RN 305 

EN;> 306 
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REAL FUNCTION KO<X) 307 

REAL 10 308 

KO ■ 0* 309 

IF(X.LT.2.> GO TO 120 310 

X2 * 2*/X 311 

KO ■ 1./SQRT<X)*EXP<-X)*<1. 25331414“. 07832358*X2 +.02189568*X2**2 312 

2 -.0l062446*X2**3+.00587872*X2**4-.0025154*X2**5+.00053208*X2**6) 313 

RETURN 314 

120 IF(X.LT.O*> RETURN 315 

T = X/3.75 316 

X2 ■ X/2. 317 

I0 C 1.+3.5156229*T**2+3.0899424*T**4+1.2067492*T**6+.2659732*T**8 318 

2 +.0360768*7**10 +.0045813*T**12 319 

KO— ALOG(X2)*I0“.57721566+.4227842*X2**2+.2306976*X2**4+ 320 

2 • 0348 859*X2**6+. 00262698*X 2**8+ .000 10 75*X 2**10+. 0000074*X2** 12 321 

RETURN 322 

END 323 

REAL FUNCTION K1(X) 324 

REAL II 325 

K1 « 0. 326 

IF(X.LT,2.) GO TO 130 327 

X2 * 2./X 328 

K1 » 1 . /SORT ( X ) *EXP ( -X )* ( 1 .25331414+. 2349861*X2”.0365562*X2**2 + 329 

2 .01504268*X2**3-.00780353*X2**4+,00325614*X2**5-.00068245*X2**6) 330 

RETURN 331 

130 IF(X.LT.O.) RETURN 332 

T = X/3.75 333 

X2 = X/2. 334 

II * X*< .5+.87890594*T**2+.51498869*T**4+.15084934*T**6 + 335 

2 .02658733*T**8+.00301532*T**10+.000324U*T**12> 336 

K1 * I 1*AL0G( X2 )+l./X*( l.+. 15443 144*X2**2“. 672 78579*X2**4“ 337 

2 • 18 156897*X 2**6“ .0191 940*X2**8“» 001 10404*X2**10- ,0O004686*X2**12 > 338 

RETURN 339 

END 340 
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subroutine KKP( X*P*CAY tCAYP ) 341 

COMMON /SUB/ PI*HALFPI 342 

DATA(KSWITCH«1) 343 

REAL K1*K0 344 

345 

346 

K SUB I NU (NU 2) SEE PAGE 5* 11(C) 347 

348 

CAY * 0*0 349 

CAYP « 0*0 350 

IF ( P *GT* 0*0000001) GO TO 200 351 

CAY = KOCX) 352 

CAYP » - K1IX) 353 

RETURN 354 

200 2 * X / P 355 

I F ( 2 .LT. 1*0) RETURN 356 

GO TO (250*350) *KSWITCH 357 

358 

C INITIALI2E AB FOR KIDKI 359 

360 

250 KSWITCH = 2 361 

CALL AB(1*DUM1*DUM2*DUM3*DUM4*DUM5*DUM6*0UM7) 362 

350 CALL KIDKI(P*2*CAY*CAYP) 363 

RETURN 364 

END 365 
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SUBROUTINE KIDKI (NU*Z*CAY*CAYP ) 366 

367 

RETURNS EXP(PI Z/2) K AND K PRIME 368 

369 

DIMENSION A(2) *B<2) *CV<2) *DV<3) 370 

DIMENSION GAM(6) *C<6) »S(6) *CP(6) 371 

COMMON /SUB/ PI.HALFPI 372 

REAL NU*NUSO*NUFTH 373 

DOUBLE ZETA*TZET*U»ZSQ 374 

DATA ( GAM = 2.6789385347, 1.3541179394,1.0*0.8929795116, 375 

* 0.9027452929,1.0)* 376 

* (S = 0.866025403784*0.866025403784,0.0,-0.866025403784, 377 

* -0.866025403784,0.0), (PI = 3.14159265359) 378 

DATA ( HALFPI = 1.570796327) 379 

CONST = SQR T ( 2 • 0 ) * PI 380 

ZSQ = Z**2 - 1.0 381 

TWOTHD = 2. 0/3.0 382 

I F ( (NU**TWOTHD * ZSQ) .LT.0.01 )GO TO 200 383 

384 

BALOGH* S FORMULA SFE PAGF 5, 11(C) 385 

386 

NUSQ = NU**2 387 

NUFTH = NUSQ** 2 388 

ONETHD = 1.0/3. 0 389 

U = DSQRT(ZSQ) 390 

I F ( Z .GT. l.OOD GO TO 10 391 

TZET = 0.0 392 

DO 50 I = 1*6 393 

11=2*1+1 394 

TZET = TZET + (-l.o')**( 1 + 1 ) * U**II/II 395 

50 CONTINUE 396 

GO TO 15 397 

10 YB = AC0S(1.0/Z) 398 

TZET = U - YB 399 

15 ZETA = (1.5 * TZET) **TWOTHD 400 

XI * NU**TWOTHD * ZETA 401 

20 CALL AB(2*NU, ZETA, ZSQ, A*B*CV,DV) 402 

CALL AIRYCXI ,AI ,AIP) 403 

ATMP * A( 2) 404 

2 I F ( ABS (AI * A ( 1 ) + A I P * B( 1 ) ) .LE. ABS (AI) ) GO TO 1 405 

A ( 1 ) = B( 1 ) * ATMP = B ( 2 ) * 0.0 406 

GO TO 5 407 

1 I F ( ABS ( AI*A ( 2 ) + AIP *B(2) ) .GT. ABS (AI * A(l) + AIP*B<1)) ) 408 

* ATMP * B ( 2 ) = 0.0 409 

5 TA = AI * (1.0 + A ( 1 ) + ATMP) 410 

6 TB = AIP * ( B ( 1 ) + B ( 2 ) ) 411 

4 I F ( ABS <AIP*DV(2) + AI * CV(1)) .LE. ABS (AlP*DV(l)) ) GO TO 3 412 

CV(1) = CV( 2 ) = DV ( 2 ) = DV ( 3 ) = A(2) = 0.0 413 

GO TO 7 414 

3 I F ( ABS ( AIP*(DV(3)-A(2) )+AI*CV(2) ).GT.ABS ( AI P*DV( 2 ) +AI*CV( 1 ) ) ) 415 

• A ( 2 ) = CV ( 2 ) = DV( 3 ) = 0.0 416 

7 TPA = AI *(CV(1) + CV( 2 ) ) 417 

8 TPB = AIP * ( DV ( 1 ) + DV( 2 ) + (DV(3) - A(2) )) 418 

TPRE = CONST 419 
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TMPRE* (ZETA/2S0) **0.25 *20 

PRE * TPRE * TMPRE / (NU**ONETHD) *21 

PREP * TPRE / (Z * <NU«*TWOTHD) * TMPRE) *22 

CAY * PRE * tTA ♦ TB) *23 

CAYP * PREP * C TPA ♦ TPB) *2* 

RETURN *25 

C *26 

C ERDELYI FORMULA SEE PAGE 7* 1 1(C) *27 

C *28 

200 EX ■ NU * <Z - 1.0) *29 

EX2 * EX**2 * 30 

EX3 * EX * EX2 *31 

EX4 * EX2**2 *32 

C(l> * 1.0 *33 

C ( 2 ) * EX *3* 

C ( 3 ) = 0*5 * EX2 + 0.05 *35 

C(4) * (EX/3.0) * (0.5 * EX2 + 0*2) *36 

C ( 5 ) = (( ( EX4 + EX2) / 6.0) ♦ 1*0 / 70.0) * 0.25 *37 

C ( 6 ) = (EX/60.0)* (0.5 * EX4 + EX2 ♦ 43.0 / 80.0) *38 

CP(1) » 0.0 *39 

CP ( 2 ) * 1*0 **° 

CP ( 3 ) = EX **1 

CP ( 4) * 0.5 * EX2 + 1./ 15. **2 

CP ( 5 ) = (EX3/6.) + (EX/12.) **3 

CP ( 6 ) » (EX4/24.) + 0.05 * EX2 + 43./4800. *** 

SUM =0.0 ** 5 

SUMP * 0.0 446 

DO 300 MM * 1.6 **7 

M = MM - 1 **8 

TS ■ ( — 1. ) **M * S(MM) * GAM ( MM) 449 

TNUZ = NU*Z/6. *50 

TX * TNUZ **(-( (M+l.)/3.) ) *51 

SUM = SUM + TS * TX * C(MM) *52 

SUMP * SUMP + TS *(CP(MM) * TX - ( (M+l. ) /18. )*C(MM)*TNUZ**(-( (M+4. 453 

1 ) /3. ) ) ) *5* 

300 CONTINUE *55 

CAY = SUM / 3.0 *56 

CAYP = SUMP / 3.0 *57 

RETURN *58 

END *59 
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SUBROUTINE AB ( LSW I TCH *NU *ZETA *ZSQ # ACC *BCC *CVC *DVC ) 

SEE PAGE 7 

DIMENSION ACC<2) *BCC(2) *CVC(2) *DVC<3) 

DOUBLE A(5)*B(5)*AC(2)*BC<2)*DV(3)*U(5)»VV(12)*D(12)*ZT(12)*VP(5)* 

1 CV(2> *NU13*NU23*DN*ZSQ*V*TEMPVP* TEMPcP* ZETA* 

2 V13 * DDt ZNU* RCOEF* COEF* PRODA* PRODB * PRODC* DVP* DVV 
GO TO ( 100*200) *LSWITCH 

INITIALIZE A»S AND B*S 

100 A(l) = 1.0 
B ( 1 ) = 1.0 
DO 150 N = 2*5 
IS * N - 1 
DN * 1.0 
ME * 2 * IS 
DO 130 M ■ 1*ME 

130 DN * DN * (2.0 * IS + 2.0 * M - 1) 

A ( N ) = DN / ( FACT (IS) * 144.0**IS) 

B(N) = -((6.0 * IS + 1.0) / (6.0 * IS - 1.0)) * A ( N ) 

150 CONTINUE 
RETURN 

LSW ITCH = 2 

200 IF ( ZSQ .GT. 0.0) GO TO 201 
PRINT 600 

600 FORMAT ( 1H0**AB ERROR ZSQ ■ **F.20.9) 

CALL EXIT 

201 V = 1.0 / DSQRT(ZSQ) 

NU13 = NU**( 1./3. ) 

NU23 = NU** ( 2 • /3 . ) 

V13 = V / NU 1 3 

DD = 39813120* 

ZNU = NU * ZETA**1 • 5 
RCOEF = ZNU**(l./3.) 

COEF = 1. /RCOEF 
PROOA = 1.0 

PRODB = 1.0 

PRODC = 1.0 

DO 300 M = 1*12 
PRODA = PRODA * V13 

PRODB = PRODB * NU23 

PRODC = PRODC * ZNU 

VV( M ) * PRODA 

D(M) = PRODB 
ZT ( M ) = PRODC 
300 CONTINUE 
U(l) = 1.0 

U( 2 ) = (1./8.)* ( ( 5«/3. ) * VV ( 3 ) + VV(1) /D<1)> 

U ( 3 ) = ( ( 9./128. ) *VV( 2 ) /D<2) +(77./l92.) * VV(4) / D(l) 


460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 
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* ♦ (385.0/1152.) * VV ( 6 ) ) 514 

U(4) « ((75./1024.) * VV(3) / D(3) ♦ ( 4563./5120.) * VV(5)/D(2) 515 

* ♦ ( 17017. /9216.) * VV( 7) / D<1) ♦ 1 05O85./82944. ) * VV(9)) 516 

U(5) • ( (4465125. /DD) * VV(4> / D(4) ♦ (94121676. /DD) *VV(6)/D(3) 517 

* + ( 349922430. 0/DD) * VV(8) / D(2) ♦ (446185740. /DD ) * VV(10) 518 

#/D(l) + ( 105910725. /DD) * VV( 12 ) ) 519 

DVP ■ 414720* 520 

VP(1) » 1.0 521 

VP(2) * -U3./8.) * VV( 1 )/D( 1 ) ♦( 7./24. ) * VV(3)) 522 

VP( 3 ) « -((15./128. )*VV( 2)/D(2) ( 33. /64. )*VV(4) / D(l) 523 

1 ♦ (455./ 1152. )*VV(6) ) 524 

VP( 4 ) * -((42525. /DVP)*VV(3)/D(3) + (451737./DVP) *VV(5)/D(2) 525 

1 +(883575. /DVP)*VV(7)/D( 1) +(475475./D VP )*VV(9)i 526 

DVV * 39813120. 527 

VP( 5) * -( (5740875 ./DVV) * VV(4) / D(4) ♦ ( 1U234708*/DVV) * 528 

1 VV( 6 ) / D(3» + (396578754. /DVV) * VV(8) / D(2) + 529 

2 (493152660. /DVV) * VV(IO) / DU> ♦ ( 202076875. /DVV) 530 

3 * VV (1 2 ) ) 531 

DV( 1 ) = 1.0 532 

DO 450 J*2*3 533 

JJ • 2 * J - 1 534 

DV(J) = 0.0 535 

DO 400 I * 1#JJ 536 

II • JJ - I + 1 537 

TEMPVP = (-l.)**(I+l) * A ( I ) * VP(II) 538 

IF(I.GT.l) TEMPVP * TEMPVP/ZT ( 1-1 ) 539 

400 DV( J) = DV( J ) + TEMPVP 540 

450 CONTINUE 541 

CV(1) = RCOEF * ( BID * VP<2) - B(2) * VP(1) / ZT(1)) 542 

CV( 2 ) * 0.0 543 

DO 475 J = 1*4 544 

JJ ■ 5 - J 545 

TEMPCP = (-l.)**(J-l) * B(J) * VP(JJ) 546 

I F ( J .GT, 1) TEMPCP = TEMPCP / ZT(J-l) 547 

475 CV( 2 ) = CV( 2 ) + TEMPCP 548 

CV( 2 ) = RCOEF * CV( 2 ) 549 

AC ( 1 ) * B(l> *U ( 3 ) - B ( 2 ) * U( 2 ) / ZT(1) + B(3) * U(l) /ZT(2) 550 

AC ( 2 ) = BID * U ( 5 ) - B ( 2 ) * U(4)/ZT(1) + 013) * U(3) / ZT(2) 551 

* - B ( 4 ) * U(2) / ZT ( 3 ) + B ( 5 ) * U(l) / ZT(4) 552 

BC(1) = COEF * ( A ( 1 ) * U(2 ) - A(2) * Ud) / ZT(1)) 553 

BC(2 ) = COEF * ( A( 1 ) * U(4) - A(2) * U(3) / ZT(1) + A(3) * U(2) 554 

* / ZT ( 2 ) - A(4) * U(l) / ZT( 3 ) ) 555 

DO 485 LU = 1*2 556 

ACC(LU) * AC(LU) 557 

BCC(LU) = BC (LU) 558 

CVC(LU) = CV(LU) 559 

DVC(LU) = DV(LU) 560 

485 CONTINUE 561 

DVC ( 3 ) « DV( 3 ) 562 

RETURN 563 

END 564 
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C 

C 

C 


SUBROUTINE AIRY (X*ERI*ERIP) 565 

566 

TABLES FROM ABRAMOWITZ AND STEGUN SEE PAGE 7 567 

568 

DIMENSION AI(101*2)»RZETA(21),FG(21*2)»XX(101) 569 

DATAUX = .00» *01* *02* .03* .04* .05, .06. .07* .08* .09* 570 

•10* .11* .12* .13* .14* .15* .16* .17* .18* .19* 571 

•20* .21* .22* .23* .24, .25* .26, .27* .28, .29* 572 

.30* .31* .32* .33* .34, .35, .36* .37, .38, .39* 573 

.40* .41* .42* .43* *44, .45* .46, .47* .48* .49* 574 

.50* .51* .52* .53* .54, .55* .56, .57* .58, .59* 575 

*60 • .61* .62* .63* .64* .65* .66* .67, .68* .69* 576 

.70. .71, .72, .73, .74, .75, .76, .77, .78, .79, 577 

. 8o » .81* .82* .83* .84* .85, .86* .87* .88* .89* 578 

.90* .91, .92, .93, .94, .95, .96, .97* .98, .99* 1.0) 579 


DATA <<AIU>* 

1=1*95) = 




580 

# 

0.35502805, 

0.35243992* 

0.34985214* 

0.34726505* 

0.34467901, 

581 

* 

0.34209435, 

0.33951139* 

0.33693047, 

0.33435191, 

0.33177603* 

582 

* 

0.32920313* 

0.32663352* 

0.32406751, 

0.32150538* 

0.31894743* 

583 

# 

0.31639395, 

0.31384521* 

0.31 130150 * 

0.30876307, 

0.30623020* 

584 

* 

0.30370315* 

0.30118218* 

0.29866753* 

0.29615945* 

0.29365818* 

585 

# 

0.29116395, 

0.28867701* 

0.28619757* 

0.28372586, 

0.28126209, 

586 

# 

0.27880648 , 

0.27635923* 

0.27392055, 

0.27149064, 

0.26906968, 

587 

* 

0.26665787, 

0.26425540 * 

0.26186243, 

0.25947916* 

0.25710574, 

588 

# 

0.25474235 , 

0.25238916, 

0.25004630* 

0.24771395* 

0.24539226, 

589 

* 

0.24308135, 

0.24078139, 

0.23849250* 

0.23621482, 

0.23394848, 

590 

# 

0.23169361 , 

0.22945031, 

0.22721872* 

0.22499894, 

0.22279109, 

591 

* 

0.22059527, 

0.21841158, 

0.21624012, 

0.21408099, 

0.21193427, 

592 

* 

0.20980006, 

0.20767844* 

0.20556948, 

0.20347327, 

0.20138987, 

593 

# 

0.19931937, 

0.19726182, 

0.19521729, 

0.19318584, 

0.19116752, 

594 

* 

0.18916240* 

0.18717052, 

0.18519192 , 

0.18322666, 

0.18127478, 

595 

# 

0.17933631* 

0.17741128, 

0.17549975, 

0.17360172, 

0.17171724, 

596 

* 

0.16984632* 

0.16798899* 

0.16614526* 

0.16431516, 

0.16249870, 

597 

# 

0.16069588, 

0.15890673* 

0.15713124* 

0.15536942, 

0.15362128, 

598 

# 

0.15188680, 

0.15016600, 

0.14845886* 

0.14676538, 

0.145Q8555) 

599 

DATA (<AI(!>» 

1=96,101) = 




600 

* 

0.14341935, 

0.14176678, 

0.14012782* 

0.13850245* 

0*13689066, 

601 

# 

0.13529242) 





602 

DATA ((AKI), 

1=102,196) = 




603 

# 

-0.25881940* 

-0.25880174, 

-0.25874909, 

-0.25866197, 

-0.25854090, 

604 

* 

-0.25838640* 

-0.25819898, 

-0.25797916, 

-0.25772745, 

-0.25744437, 

605 

4f 

-0.2571304 2 » 

-0.25678613, 

-0.25641200* 

-0.25600854, 

-0.25557625, 

606 

* 

-0.25511565, 

-0.25462724, 

-0.25411151, 

-0.25356898, 

-0.25300013, 

607 

* 

-0.25240547, 

-0.25178548* 

-0.25114067* 

-0.25047151, 

-0.24977850, 

608 

* 

-0.24906211, 

-0.24832284* 

-0.24756115* 

-0.24677753* 

-0.24597244, 

609 

* 

-0.24514636, 

-0.24429976, 

-0.24343309* 

-0.24254682* 

-0.24164140, 

610 

# 

-0.24071730* 

-0.23977495, 

-0.23881481, 

-0.23783731* 

-0.23684291, 

611 

* 

-0.23583203* 

-0.23480512* 

-0.23376259, 

-0.23270487, 

-0.23163239, 

612 

* 

-0.23054556* 

-0.22944479* 

-0.22833050* 

-0.22720310* 

-0.22606297, 

613 

* 

-0.22491053, 

-0.22374617, 

-0.22257027, 

-0.22138322, 

-0.22018541, 

614 

* 

-0.21897720* 

-0.21775898, 

-0.21653112, 

-0.21529397* 

-0.21404790, 

615 

* 

-0.21279326, 

-0.21153041, 

-0.21025970* 

-0.20893146* 

-0.20769605, 

616 

* 

-0.2Q640378, 

-0.20510500* 

-0.20380004* 

-0.20248920* 

-0.20117281, 

617 

* 

-0.19985119, 

-0.19852464, 

-0.19719347, 

-0.19585798, 

-0.19451846, 

618 
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* -0# 1931752 1 * -0* 19182851 » -0*19047865# -0*18912591* -0*18777055# 619 

* -0.18641286* -0.18505310# -0*18369153* -0*18232840# -0.18096398, 620 

* -0*17959851* -0.17823223* -0.17,^6539* -0.17549823* -0,17413097, 621 

* -0,17276384# -0,17139708* -0,1700309ft* -0,16866551* -0,16730113) 622 

OATA ( ( AI ( I ) * I * 197*202) * 823 

* -0,16593797, -0,16457623* -0*16321611* -0,16185781, -0,16050153* 624 

* -0*15914744) 625 

DATA ( RZETA * 1,5* 1.4, 1,3#1,2, 1.1, 1.0* 0.9* 0,8, 0.7, 0.6# 626 

* 0,5, 0*45, 0,4* 0,35, 0*3* 0,25* 0,2# 0,15, 0.1# 627 

* 0*05* 0*0) 528 

DATA (FG * 0.527027, 0.528783* 0.530601* 0.532488,0.534448, 629 

* 0.536489, 0.538618* 0.540844, 0.543180, 0.545636, 630 

* 0.548230, 0.549584, 0.550980, 0.552421, 0.553912, 631 

* 0.555456, 0.557058*0.558724# 0.560462, 0.562280* 632 


* 

0.564190* 





* 

0.619954, 

0.617156# 

0.614275, 

0.611305, 

0.608239, 

* 

0.605068, 

0.601782, 

0.598372* 

0.594823, 

0.591120* 

* 

0.587245 , 

0.585235* 

0.583174, 

0.581056, 

0.578878, 

# 

* 

0.576635, 

0.564190) 

0.574320* 

0.571927, 

0.569448, 

0.566873, 


NXX = 101 
NRZ =21 

I F ( X .GE. 0*0) GO TO 100 
PRINT 600* X 

600 FORMAT ( 1H0**AIRY X = *,E2o*9) 

CALL EXIT 

100 I F ( X .GT. 1.0) GO TO 300 
DO 110 J = 1 »NXX 
I F ( X - XX( ) 115,105,110 
105 ERI = A I ( J , 1 ) 

ERIP = AI ( J*2) 

RETURN 

110 CONTINUE 

115 JSTART = J - 2 

I F ( JST ART .EQ. 0> JSTART » 1 

IF( JSTART .EQ. (NXX-2)) JSTART = NXX - 3 

ERI = 0.0 

ERIP = 0.0 

DO 130 K = 1*4 

KK « K + JSTART - 1 

DN = 1.0 

DD = 1.0 

DO 120 J = 1*4 

I F ( J .EQ. K ) GO TO 120 

JJ = J + JSTART - 1 

DN = DN * (X - XX ( JJ) ) 

DD = DD * <XX(KK) - XX (JJ)) 

120 CONTINUE 

TEMP = DN/DD 

ERI = ERI + TEMP * AI(KK*1) 

ERIP = ERIP + TEMP * AI(KK,2) 

130 CONTINUE 
RETURN 


633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 
667 
666 

669 

670 

671 

672 


C 


n n 
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X *GT. 1*0 673 

674 

300 ZETA * (2. 0/3*0) * X**1.5 675 

RZET ■ 1.0/ZETA 676 

XFOURTH ■ X**0.25 677 

XPONE - 0*5 * EXP ( -ZETA ) 678 

DO 310 M * l.NRZ 679 

IF(RZET - RZETA(M)) 310*305. 315 687 

305 F « FG(M.l) 681 

G ■ FG(M*2) 682 

GO TO 350 683 

310 CONTINUE 684 

315 JSTART * M- 2 685 

I F ( JST ART *EQ. 0> JSTART » 1 686 

IF( JSTART • EO* ( NRZ - 2)) JSTART * NRZ - 3 687 

F = 0*0 688 

G = 0.0 689 

DO 330 K * 1 »4 690 

KK * K + JSTART - 1 691 

DN * 1*0 692 

DD * 1.0 693 

DO 320 J = 1 *4 694 

IF ( J .EQ. O GO TO 320 695 

JJ = J + JSTART ~ 1 696 

DN = DN * (RZET - RZETA(JJ)) 697 

DD * DD * (RZETA(KK) - RZETA(JJ)) 698 

320 CONTINUE 699 

TEMP * DN / DD 700 

F = F + TEMP * FG( KK * 1 ) 701 

G = G + TEMP * FG ( KK *2 ) 702 

330 CONTINUE 703 

350 ERI = XPONE * F / XFOURTH 704 

ERIP « (-XPONE) * G * XFOURTH 705 

RETURN 706 

END 707 


FUNCTION FACT ( N ) 708 

TYPE DOUBLE FN 709 

IF ( N ) 100*101*102 710 

100 PRINT 300 *N 711 

300 FORMAT ( 1H0* 10X »*NEGAT IVE FACTORIAL ARGUMENT** 15) 712 

CALI EXIT 713 

101 FACT * 1.0 714 

RETURN 715 

102 FN * 1.0 716 

DO 105 I * 1 *N 717 

105 FN * FN * I 718 

FACT = FN 719 

RETURN 720 

END 721 
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B(£,e) is calculated by a FORTRAN sub-program 
SUBROUTINE BFINAL. 
a typical call to BFINAL would be: 



40 

41 

42 

42 

43 
43 


i 


ten. 


oncD>^a>-Nj^ 4 >^wfOM 03 i\) h 
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SUBROUTINE BFINAL(P*EPS*B> 

COMMON /TABLE/ EPSILON ( 15 *3 ) *PF( 65 * 3 ) *BF( 65. 15 *3 ) *NEPS< 3 ) »NeP ( 3 ) . 

* PLASTO) 

DATA(NEPS=9*!2*8) ,(NEP«65*43*38) *(PLAST«12.5*30.0*5.0) 

DATA! EPSILON ■ 1.0001* l.OOl* 1.01* 1.03* 1.05* 1.1* 1.2* 1.25* 

1.3* 6(0.0)* 1.3* 1.35* 1.4* 1.45* 1.5* 1.55* 1.6* 

1.65* 1.7* 1.8* 1.9* 2.0* 3(0.0)* 2.0* 3.0* 5.0* 

8.0* 10.0* 20.0* 50.0* 100.0*7(0.0) ) 

DATA ( PF = 0*0*0.02* 0.04* 0.06* 0.08* 0.1* .12* .14. .16* .18* .2* 
.4* .6* .8* 1.0* 1.2* 1.4* 1.6* 1.8* 2.0* 2.2* 2.4* 2.6* 
2.8* 3.0* 3.2* 3.4* 3.6* 3.8' 4.0* 4.2* 4.4* 4.6* 4.8* 

5.0* 5.2* 5.4* 5.6* 5.8* 6.0* 6.2* 6.4* 6.6* b.8* 7.0* 

7.2* 7.4* 7.6* 7.8* 8.0* 8.2* 8.4, 8.6, 8.8* 9.0* 9.2* 

9.4* 9.6, 9.8* lO.O • 10.5* 11.0* 11.5* 12.0* 12.5* 

0*0* .02* .04* .06 * .08* .1* .2* .3* .4* .5* .6* .7* .8* 
.9* 1.0* 1.5* 2.0* 2.5* 3.0* 3.5* 4.0* 4.5* 5.0* 5.5* 
6*0* 6.5* 7.0* 7.5* 8.0* 8*5* 9.0* 9.5* lO.O* 12.0* 14.* 
16.0* 18.0* 20.0* 22.0* 24.0* 26.0* 28.0* 30.0* 22(0.0)* 
.0025* .005* .0075* .01* .02* .04* .06* . 08 * .1* .12* 
.14* .16* .18* .2* .4* .6* .8* 1.0, 1.2* 1.4* 1.6* 1.8* 
2.0* 2.2, 2.4* 2.6, 2.8, 3.0* 3.2, 3.4* 3.6* 3.8* 4.0* 

4.2* 4.4* 4.6* 4.8* 5.0* 27(0.0) ) 

C 

C SMOOTHED B-TABLE 1 

C 1.0001 LE EPSILON LT 1.3 SEE PAGE 10 (A) 

C 


DATA ( <BF( I )» 1=1,130)= 


* 

0.00000*- 

.01705, 

-.03410,- 

.05115, -.06820* 

-.07641*- 

•08180* 

-.08719* 

* 

-.09258,- 

• 09643 * 

-.09885,- 

.10620, -.09039* 

-.06446,- 

.03479* 

-.0044j, 

* 

.02487, 

.05211 , 

.07669, 

.09842, ,11734, 

. 13292, 

.14572* 

.15396* 

* 

.16220* 

.17044, 

.17847, 

.18001* .18155* 

.18309, 

•18462, 

.18335* 

* 

.18132* 

.17929, 

.17726, 

.17439, .17082* 

• 16725* 

•16368, 

.15999, 

* 

.15602, 

.15205 , 

.14808* 

.14412, .14032, 

.13652, 

.13272* 

•12892, 

# 

.12545, 

• 1 2 2 06 » 

.11867, 

.11528, .11217, 

.10925, 

•10634, 

.10342*, 

* 

• 10066* 

.09819* 

.09573, 

.09327, ,08775* 

.08259, 

.07821* 

.07406* 

* 

.07042* 







* 

o.ooooo*- 

•01658* 

-.03316*- 

,04974, -.06633 , 

-.07678,- 

.08187, 

-.08696* 

* 

-.09204*- 

.09694, 

-.09906,- 

. 10605 *-.09016 , 

-.06439,- 

•03568* 

-.00743, 

* 

.02082, 

.04907, 

.07658, 

•09368, .11079, 

.12789* 

.14500* 

.15587, 

* 

•16258* 

.16930* 

.17602* 

.18190* .18234* 

.18279* 

.18323* 

•18367, 

* 

• 18206 » 

.17941, 

.17677, 

.17412, .17125, 

. 16744, 

.16363, 

.15983, 

* 

.15602* 

.15211, 

.14816, 

.14420, .14025* 

.13639, 

.13275, 

.12912, 

# 

.12548* 

.12185, 

.11860, 

.11545, .11229, 

. 10914, 

.10615, 

.10349, 

* 

.10084 » 

.09818, 

.09552, 

.09325, .08771, 

.08271, 

.07812, 

.07413* 

* 

•07031 ) 







DATA ( <3F( I ) 

,1 = 131, 

260) = 





# 

0.00000*- 

.03204, 

-.04764,- 

.06048, -.07022, 

— . 07756 , — 

.08370* 

-.08895* 

* 

-.09334*- 

.09703, 

-.10009,- 

. 10132 , — .08867 »■ 

-.06247,- 

.03269, 

-.00235, 

* 

.02751, 

• 05445 * 

.07896, 

.10108, .12001, 

.13571, 

.14919, 

•16052* 

* 

. 16883 * 

.17542, 

• 18050 * 

.18387, .18567, 

• 18652, 

.18647, 

.18530, 


.18353, 

.18126, 

.17847, 

.17525, .17179, 

.16813, 

. 16425 , 

.16028, 

* 

. 1 5626 * 

.15221, 

.14816, 

.14415, .14018, 

.13630* 

.13250* 

• 12878* 

# 

.12517, 

.12169, 

.11830* 

.11501* .11187, 

• 10884* 

.10590* 

• 10308* 

# 

. 10040 * 

.09780* 

.09530, 

.09292, .08738, 

.08239, 

.07789, 

.07383* 


001 

002 

003 

004 

005 

006 

007 

008 

009 

010 

011 

012 

013 

014 

015 

016 

017 

018 

019 

020 
021 
022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 

038 

039 

040 

041 

042 

043 

044 

045 

046 

047 

048 

049 

050 

051 

052 

053 

054 


* * 
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# 

•07018* 




# 

0.00000*- 

-.03243* 

-.04897* 

-.06121, -.07094, -.07839, -.08438, -.08948 

* 

-.09378,- 

.09745, 

-.10046, 

-.10559, -.08752, -.06003, -.02889, .00259 

# 

.03279, 

•06048* 

•08540* 

.10734, .12604, .14205, .15507* .16571 

# 

.17363* 

•18006, 

.18491, 

.18755, .18914, .18973, .18913* .18769 

# 

•18569* 

.18312, 

.17998* 

.17657, .17290, .16898, .16493, .16081 

# 

•15663* 

.15245, 

.14829, 

.14416, .14013, .13617, .13229, .12853 

* 

.12489, 

•12134, 

.11791, 

.11464, .11146, .10839, .10547, .10266 

* 

.09995* 

.09736, 

.09489, 

.09250, .08696* .08201, .07755, .07353 

* 

.06989 ) 




DATA ( ( BF( I ) 

*1=261, 

390) = 


# 

0.00000*- 

.03243, 

-.04886, 

-.06 139, -.07 10 5, -.07863, -.08475, -.08987 

# 

-.09428*- 

.09783, 

-.10062* 

104 59, -.08565 *-.05695 *-.02484, .00741 

* 

.0379 1 , 

•06613 , 

.09122, 

.11257, .13122, .14754, .16003, .17029 

* 

.17874, 

.18444, 

.18853* 

.19139, .19248, .19251, .19178, .19001 

# 

.18762* 

.18479, 

.18143* 

.17775, .17386* .16975, .16553, .16125 

* 

.15694* 

.15264, 

.14837* 

.14422, .14007, .13601, .13216, .12832 

# 

.12458, 

. 12110* 

.11762* 

.11427, .11116, .10806, .10508, .10234 

# 

.09959* 

.09697, 

.09456 * 

.09214, .08666, .08172, .07727, .07325 

* 

•06962* 




# 

o.ooooo*- 

.03198, 

-.04914, 

-.061 56, -.071 39 ,-.0782 7, -.08509, -.08984 

# 

-.09452*- 

.09758, 

-.10057, 

-. 10242 ,-.08039, -.04901, -.01471, .01921 

# 

.05116, 

.08014, 

.10541, 

.12730, ,14588, .16105, .17213, .18282 

# 

• 18908, 

•19500, 

,19757* 

.19986, .19973. .19938, .19736, .19517 

* 

.19188* 

. 1 8846 , 

.18419* 

.17992, .17564, .17132, .16669, ,16206 

* 

.15743, 

.15284, 

.14847, 

.14410, .13972, .13545, .13158, .12771. 

* 

.12384* 

.12010* 

.11679, 

.11347, .11016, .10699, .10420* ,10141i 

# 

.09862* 

.09596, 

.09363* 

.09130, .08578* .08090* .07659, .07257 ( 

* 

.06910 ) 





data 
* 0 


* 

* 

* 

# 

* 

# 

* 


# 

# 

# 

* 

# 

* 

* 

* 

DATA 

* 0 
* - 

* 

* 

# 


( <BF<I> 
00000 *-* 
09397*-. 
07687* . 
20880 * • 
19859, • 
15764* . 
12228* • 
09692, . 
06816* 
ooooo*-* 

09336,-. 
08909, . 
21648, . 
20093* • 
15731, . 
12106, • 
09591, • 
06782 


) 


(BF( I ) 

ooooo,-* 

09239,-. 
lQlOl* • 
22618, • 
2 q290* • 


*1=391 

03157 
09702 
10687 
21265 
19388 
15285 
11833 
09447 

03161 

09613 

11954 

22022 

19581 

15207 

11729 

09354 

*1= 52 

03158 
09511 
13153 
22766 
19742 


,520)= 

-.04887 

-.09933 

.13262 

.21368 

•18886 

•14806 

.11476 

.09202 

-.04785 

-.09824 

.14520 

.22016 

.19035 

.14712 

.11391 

.09139 

1* 975) 
-.04792 
-.09714 
.15722 
.22726 
.19147 


-.06151 

-.09600 

.15420 

.21394 

.18368 

.14327 

.11163 

.08956 

-.06105 

-.09242 

.16656 

.21914 

.18474 

.14234 

.11091 

.08923 


-.06047 

-.08722 

.17827 

.22420 

.18553 


-.07131 

-.06876 

.17012 

.21234 

.17842 

.13848 

•10850 

.08438 

-.07111 

-.06211 

.18363 

.21786 

.17914 

.137-^7 

.10791 

.08384 


-.07027 

-.05547 

.19501 

.22114 

.17958 


-.07893,- 

-.03258, 

.18540, 

.21015, 

.17311, 

•13413, 

.10537, 

.07954, 

-.07845,- 

-.02415, 

.19698, 

.21410, 

.17357, 

.13325, 

.10491, 

.07902, 


-.08503,- 

.00592* 

.19524, 

. 20680 , 

.16784, 

•13018, 

•10224, 

.07539, 

.08454,- 

.01600* 

.20691, 

.21035, 

.16801, 

.12899, 

.10191, 

.07506, 


-.09004 

.04296 

.20403 

•20302 

•16261 

•12623 

.09938 

.07153 

-.08952 

•05450 

.21274 

.20605 

.16254 

.12482 

.09891 

.07111 


-.07831,-. 08394 ,-.088‘ 
.01525, .02599, .065' 
.20746, .21648, .222! 
.21808, .21354, .2o8J 
.17382, .16808, .162: 


055 

056 

057 

058 

059 

060 
061 
062 
063 

064 

065 

066 
067 
066 

069 

070 

071 

072 

073 

074 

075 

076 

077 

078 

079 

080 
081 
082 

083 

084 

085 

086 

087 

088 

089 

090 

091 

092 

093 

094 

095 

096 

097 

098 

099 

100 
101 
102 

103 

104 

105 

106 

107 

108 


nnnn 
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# 

•15689. • 

15173, .14657, .14147, .13701, .13254 

, .12808 

, • 12408 , 

109 

# 

• 12028 . o 

11648, .11286, .10966, .10646, .10326 

, .10052 

♦ .09782# 

110 

* 

•09513. • 

09265, .09038, .08811, .08332, ,07899 

, .07467 

, .07068, 

111 

# 

•06766. 390(0.0) ) 






112 









113 

ACTUAL B-TABLE 2 






114 


1.3 LE EPSILON LE 2.0 



SEE PAGE 11, (B) 

115 









116 

DATA ( < BF ( I > 

*1 = 976,1105)= 





117 

* 

o.ooooo* 

-.16348, 

-.26064* 

-.33995, 

- .40848 , 

-.46904, 

-.69532, 

118 

# 

- . 840 33 ♦ 

-.93753, 

-1.00021*' 

-1.03685,' 

-1.05300,- 

1.05263, 

-1.03870, 

119 

* 

-1.01358, 

-.77399, 

-.43202 , 

— . 0560 6 , 

.31451, 

. 6 7404 , 

1.00020, 

120 

* 

1.28943, 

1.54059, 

1.75435, 

1.93254, 

2.07771, 

2.19280, 

2.28089, 

121 

* 

2.34504, 

2.38813, 

2.41310* 

2.42231 , 

2.41830, 

2.30789, 

2.11598, 

122 

* 

1.90148, 

1.69468, 

1.50860* 

1.34756, 

1.21059, 

1.09511, 

.99792, 

123 

* 

.91601 

, 22(0.0) 

, 





124 

* 

0.00000* 

-.16105, 

-.25564, 

-.33189, 

-.39700, 

-.45399, 

-.65836, 

125 

* 

-.77968, 

-.85118, 

-.88732 » 

-.89736, 

-.88708, 

-.86093, 

-.82205, 

126 

* 

-.77294, 

-.43221 , 

-.02425 , 

.38313, 

.75521, 

1.09165, 

1.37352, 

127 

* 

1 .6044 1 , 

1.78758, 

1.92759, 

2.02956, 

2.09862, 

2.13975, 

2.15749, 

128 

* 

2.15594, 

2.13895, 

2.10943 ♦ 

2.07035, 

2.02387, 

1.79949, 

1.56531, 

129 

* 

1.35512, 

1.17858, 

1,03445 , 

.91808, 

.82377, 

.74682, 

.68343, 

130 

* 

.63029 

,22(0.0) 

) 





131 

DATA ( <BF(I) 

,1=1106*1235)= 





132 

* 

0.00000* 

-.15872, 

-.25074, 

-.32398, 

-.38579, 

-.43905, 

-.62235, 

133 

* 

-.72081, 

-.76764* 

-.77902, 

-.76442, 

-.73053, 

-.68170, 

-.62139, 

134 

* 

-.55225, 

-.13486, 

.30988, 

.71833, 

1.06787, 

1.35450, 

1.57566, 

135 

* 

1.73830* 

1.85000, 

1.91871 , 

1.95267, 

1.95875, 

1.94304, 

1.91075, 

136 

* 

1.86613* 

1.81274, 

1.75348 , 

1.69063 , 

1.62601, 

1.37300, 

1.15558, 

137 

* 

.98234, 

.84786, 

.74367, 

.66203, 

.59696, 

.54407, 

.50027, 

138 

* 

.46338, 

22(0.0) 

, 





139 

* 

0.00000* 

-.15637, 

-.24594 , 

-.31621 , 

-.37459, 

-.42437, 

-.58736, 

140 

* 

-.66323 » 

-.68678, 

-.67518, 

-.63837, 

-.58337, 

-.51489* 

-.43656, 

141 

* 

-.35112* 

.12068, 

.57783, 

.96502, 

1.27087, 

1.49834, 

1.65372, 

142 

* 

1.74935, 

1.79706, 

1.80753, 

1.78994, 

1.75194, 

1.69966, 

1.63804, 

143 

* 

1.57098, 

1.50112, 

1.43085, 

1.36155, 

1.29436, 

1.05739* 

.87589, 

144 

* 

.74131, 

.64100, 

.56498 , 

.50562, 

•45816, 

.41934* 

.38691, 

145 

* 

.35929, 

22(0.0) 

) 





146 

DATA ( (BF(I> 

,1=1236,1365)= 





147 

* 

0.00000* 

-.15414, 

-.24122 , 

-.30852, 

-.36370, 

-.41001, 

-.55316, 

148 

* 

-.60728* 

— . 60 886 , 

-.57583, 

-.51910* 

-•44552, 

-.36032, 

-•26719, 

149 

* 

-.16894, 

.33761, 

.78713, 

1.13708, 

1.38973, 

1.55437, 

1.64653, 

150 

* 

1.68219, 

1.67687, 

1.64256, 

1.58875, 

1.52276, 

1.45006, 

1.37468, 

151 

* 

1.29939, 

1.22611, 

1.15608, 

1.09002, 

1.02830* 

.32494, 

• 68056, 

152 

* 

.57768, 

.50230, 

.44509, 

•40021, 

.36396, 

.33401* 

.30878, 

153 

* 

.28720* 

22(0.0) * 






154 

* 

0.00000* 

-.15203, 

-.23663, 

-.30105 , 

-.35304, 

-.39576, 

-.51961, 

155 

* 

-.55285, 

-.53376, 

-.48120, 

-.40641 , 

-.31683, 

-.21768, 

-.11275, 

156 

* 

-.00490 » 

.51924, 

.94610, 

1.24818, 

1.44370, 

1.54793* 

1.58409, 

157 

* 

1.57217, 

1.52811, 

1.46400, 

1.38877, 

1.30869, 

1.22810, 

1.14985, 

158 

* 

1.07548, 

1.00610, 

.94206, 

•88336, 

.82993, 

.66149, 

.54718, 

159 

# 

.46707* 

.40830, 

.36335, 

.32770, 

•29868, 

.27461, 

.25415, 

160 

* 

.23659, 

22(0.0) 

) 





161 


DATA ( ( BF ( I ) 

,1=1366,1495)= 





162 


nn n n 
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# 

0.00000* 

-.14985* -•23202. 

-.29354, 

-.34236, 

-.38089, - 

•48618* 

163 

# 

-•49960* 

-.46155* -.39108* 

-.30034 * 

-.19706* 

-.08654, 

.02738* 

164 

# 

•14194* 

•66908* 1.06206* 

1.30979, 

1.44976, 

1.49976* 1 

.48929* 

165 

* 

1,44050* 

1,36941, 1,28717* 

1.20130* 

1.11 666 * 

1.03615* 

.96125* 

166 

* 

•89271* 

.83058* .77469* 

.72453* 

.67969, 

•54224, 

•45103* 

167 

« 

•38708* 

.33983* ,30332* 

•27415* 

.25023* 

.23027* 

•21329* 

168 

# 

•19868* 

22(0.0) » 





169 

# 

0*00000* 

-.14775* -.22759* 

-.28618, 

-.33177* 

-.36797, - 

.45441* 

170 

# 

-•44853* 

-.39201* -.30540* 

-.20070, 

-.08593, 

.03353* 

.15393* 

171 

# 

•27257* 

.79051. 1.14185* 

1.33750* 

1.42157, 

1.42433* 1 

.37738* 

172 

# 

1.30266* 

1.21504* 1.12405* 

1.03549, 

.95258* 

.87680* 

.80851* 

173 

* 

.74772* 

.69383* .64613* 

.60401* 

•56673, 

.45424* 

.38006* 

174 

# 

.32761* 

.28853* .25802* 

•23350* 

.21340* 

•19648* 

.18208* 

175 

* 

•16967* 

22(0.0) ) 





176 

DATA ( ( BF ( I ) 

» I * 1496 * 1625 ) * 





177 

* 

0*00000* 

-.14577, -.22311* 

-.27896, 

-.32161* 

-.35447* - 

.42275* 

178 

# 

-.39868* 

-.32531* -.22399* 

—.10728* 

•01687* 

•14306* 

•26762* 

179 

# 

.38807* 

.88694* 1.19173* 

1.33567* 

1.37008* 

1.33439* 1 

•26030* 

180 

* 

1.16880* 

1.07279* .97979* 

.89363* 

•81592* 

.74708* 

.68664* 

181 

* 

.63375* 

.58763* .54735* 

.51204* 

.48098, 

.38780* 

•32610* 

182 

* 

.28205* 

•24894* .22292* 

•20195* 

. 18466 » 

•17010* 

.15769, 

183 

# 

.14699* 

22(0.0) * 





184 

# 

0.00000* 

-.14175* -.21441* 

-.26461* 

-.30153, 

-.32711* - 

.36101* 

185 

* 

“.30308* 

-.20005* -.07401* 

.06165* 

.19890, 

.33251, 

.45958, 

186 

# 

• 57784* 

1.01654* 1.22277, 

1.27257, 

1.22903, 

1.14049, 1 

.03587, 

187 

* 

.93134* 

.83474* .74917* 

.67514* 

•61181* 

.55795, 

•51212* 

188 

# 

. 47296. 

.43933, .41027, 

.38493* 

.36271* 

.29564, 

.25033, 

189 

* 

.21746* 

.19240* .17258, 

.15653* 

• 14323* 

.13203, 

•12246* 

190 

* 

.11422* 

22(0.0) ) 





191 

DATA ( (BF(I) 

,1=1626.1950)= 





192 

* 

0.00000* 

-.13797* -.20598* 

-.25079* 

-.28178, 

-.30185* - 

•30256* 

193 

* 

-.21391* 

-.08514, .05992, 

•20816* 

.35189, 

.48656, 

.60968* 

194 

* 

.71999* 

1.08124, 1.19333, 

1.16509* 

1.07114* 

.95643* 

.84384* 

195 

* 

.74300* 

.65677, .58461* 

.52475* 

.47511* 

.43381* 

.39913* 

196 

* 

.36971* 

.34451, .32270* 

.30365* 

•28683, 

.23551* 

.20019* 

197 

* 

.17432. 

•15442, .13867* 

.12585, 

.11521, 

•10623* 

.09856* 

198 

# 

.09195* 

22(0.0) . 





199 

* 

0.00000* 

-.13667, -.19753, 

-.23765, 

-.26255, 

-.27642* - 

.24616* 

200 


-.13011* 

.01966, .17869, 

. 3340&* 

.47877, 

.60917, 

.72370* 

201 

* 

.82198* 

1.09767, 1.12793, 

1.04340* 

.92091, 

.79832* 

•69036* 

202 

* 

.60048 » 

.52749, .46863* 

.42100* 

•38209* 

.34994* 

•32298* 

203 

* 

•30010* 

.28040* .26328* 

•24824* 

• 23488* 

.19370* 

•16506* 

204 

* 

.14391* 

.12760* .11465, 

.10409* 

.09532* 

.08793, 

•08160* 

205 

# 

.07612* 

217(0.0) ) 





206 








207 

SMOOTHED B-TABLE 3 





208 


2.0 6T EPSILON LE INFINITY 


SEE PAGE 11* (Cl 


209 








210 

DATA ( C BF II> 

.1=1951*2080)= 





211 

# 

-.97421* 

-.68852* -.58707* 

-.49707 

» -.28169 

, -.12075 

• 

212 

# 

-•05982* 

-.02569, -.00017* 

.02081 

, .04057 

, .05951 

* 

213 

* 

•07828* 

•09760* .30729, 

•48640 

* .57493 

, .58957 

» 

214 

* 

.57587* 

.55837, .54506, 

.53370 

» .52582 

* .52010 

» 

215 

* 

•51617* 

•51324, .51141* 

.50958 

» .50815 

, .50714 


216 


38- 


# 

. 5o613 » 

.50547, 

.50486* 

.50430* 

.50387, 

.50344, 

217 

# 

•50312* 

•50283* 

27(0.0) * 




218 

* 

-2*61045* 

-1.37569. 

-1.02995, 

-.87639, 

-.57477, 

-.30191, 

219 

* 

-•15183* 

-.07734, 

-.02930* 

.00523* 

•03119* 

.05285, 

220 

* 

•07121* 

•08632, 

.20854, 

.31337, 

.39987, 

.46999, 

221 

* 

.50859* 

•54720* 

•55318* 

.55387, 

.55145, 

•54568, 

222 

* 

•53981, 

.53315, 

.52648, 

.52236, 

.51922, 

.51648, 

223 

# 

•51450* 

•51252* 

•51112* 

.50975, 

•50864, 

.50773, 

224 

« 

•50688* 

•50626* 

27(0.0) ) 




225 

DATA ( (BF( 

I > *1=2081 , 

2210)- 




226 

* 

-7.25985* 

-3.45485, 

-1.64059, 

-1.34423* 

-.66897, 

-.34955, 

227 

* 

-.18538* 

— .06638 » 

.00132, 

• 04141 * 

.07376, 

•10610, 

228 

* 

.12558, 

• 14254 , 

.24915, 

.31181, 

.36725, 

.41363, 

229 

♦ 

•44855* 

• 47569, 

.50559, 

.51503, 

•52354, 

.53369, 

230 

* 

.53025, 

.53039, 

.52924, 

.52442, 

•52331, 

.52033, 

231 

* 

.51719* 

.51574, 

.51319, 

.51151* 

•51033* 

.50871, 

232 

* 

.50788* 

.50703 , 

27(0.0) * 




233 

*. 

>10.39392* 

-4.19995, 

-2.58240* 

-1.95097, 

-.72320* 

-.26754, 

234 

* 

-.09334* 

.02079, 

.08494* 

.13020* 

•16444* 

.17916, 

235 

«* 

.19375* 

•20833* 

.30709* 

.35538* 

•39234* 

•42172, 

236 

* 

.44781, 

.46999, 

.47925, 

•48858, 

.50351, 

.51236, 

237 

* 

,5o9o8» 

.51178, 

.52183, 

•51642, 

.51239, 

.51401, 

238 

* 

•51480* 

•51039, 

.50923, 

.51144, 

.50871, 

•50634, 

239 

* 

.50653* 

.50672, 

27(0.0) ) 




240 

DATA ( ( BF ( 

I) *1=2211* 

2340) = 




241 

*• 

-12.2234 2, 

-4.75181, 

-2.27086, 

-1.77188, 

-.63507, 

-.21446, 

242 

* 

— .04165 * 

•06149, 

.10540 * 

.14931, 

•19322, 

•22335, 

243 

* 

.23512* 

•24688* 

.33362, 

.37564, 

.40544, 

•43029, 

244 

* 

.45182, 

.46992, 

.48265, 

.49159, 

.50053, 

.50947, 

245 

* 

.51348, 

.51474, 

.51601* 

.51728, 

.51635, 

.51507, 

246 

* 

.51378, 

.51250, 

.51134, 

.51018, 

.50902, 

.50795, 

247 

* 

•50712* 

.50629, 

27(0.0) * 




248 

*-17.18560* 

-4.93449, 

-2.51731, 

- 1 .70006 , 

-.54474, 

— •Ofti'IS , 

249 

* 

.10521* 

.21329, 

.27134, 

.29964, 

.31983, 

.3686, 

250 

* 

. 349o8 » 

.35471, 

.40376, 

.42789, 

.44491 , 

.* *837, 

251 

* 

.46939, 

.47894, 

.48697, 

.49128, 

.48934, 

.•8967, 

252 

* 

.49440 » 

• 50561 , 

•50019* 

.49603, 

.49688, 

.50334, 

253 

* 

• 5o5o9 * 

•50011* 

.49846* 

• 50036 * 

•50481, 

•50104, 

254 

* 

.49933* 

.49990* 

27(0.0 ) 

) 



255 

DATA ( ( BF( 

I ) *1 = 2341. 

2925)= 




256 

*-11.38347* 

-4.45640 * 

-2.80756, 

-1.40802, 

-.36357, 

.05302, 

257 

* 

•22618, 

.32724, 

.39072, 

•40166, 

.41157, 

.41995, 

258 

* 

.42834, 

.43319, 

.45263, 

.45965* 

.47071, 

.48885 , 

259 

* 

.49584* 

.49502, 

.49897, 

.50159, 

•49886, 

.49355, 

260 

* 

•49218, 

.49272, 

.49327, 

.49381, 

.49436, 

.49490* 

261 

* 

.49545* 

.49600* 

.49654, 

.49707, 

.49747, 

.49787* 

262 

* 

.49827, 

.49867, 

27(0.0) * 




263 

* 

-6.06082* 

-3.14788, 

-1.76772, 

-1.14861, 

—.25214, 

.12786, 

264 

* 

.26495, 

.35971, 

• 44275, 

•44560* 

•45094, 

.45623, 

265 

* 

.46195, 

•46602* 

.47481, 

.47974, 

•48212, 

•48230, 

266 

* 

•47861, 

.47598, 

.48927, 

.50561, 

•51841, 

.52185, 

267 

* 

•52032* 

•51590* 

.51700* 

.51391* 

•50035* 

.49271, 

268 

* 

.49458, 

•50436, 

.50658, 

.50119, 

•49832, 

.49866, 

269 

# 

•50321* 

•50407 , 

482(0.0) 

) 



270 




t 
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C 271 

C KK ■ TABLE NUMBER 272 

C 273 

115 KK ■ 1 274 

IF<EPS .IT. 1.3) 60 TO 120 275 

KK • 2 276 

IFtEPS.GT. 2.0) KK ■ 3 277 

120 IF(KK .NE. 21 CALL CONV(EPStP) 278 

B • BOFEP<EPS*P*PTEMP»KK) 279 

RETURN 280 

END 281 


1 


n n nn 
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FUNCTION BOFEP ( E *P .PTEMP *K ) 282 

283 

INTERPOLATION IN THREE B-TABLES 284 

K * TABLE NUMBER 285 

286 

COMMON/TABLE/ EPSILON! 15.3). PF( 65 *3). BF (65.15*3) *NEPS( 3 ) *NEP( 3 ) • 287 

* PLAST ( 3 ) 288 

DIMENSION ET ( 3 ) 289 

I DUAL * 1 290 

PTEMP * P 291 

IF(K.EQ. 2) 60 TO 5 292 

CALL K1K3 ( E*P *DUM1 .PTEMP *DUM2 * 1 *K ) 293 

5 IFtPTEMP .LE. PLAST(fC) ) GO TO 10 294 

BOFEP = BASYM(E.P.K) 295 

RETURN 296 

10 IF(K.NE.3> GO TO 20 297 

I F ( PTEMP .GE. 0.0025) GO TO 15 298 

BOFEP * 0.0 299 

RETURN 300 

15 I F ( E .LE. 100.0) GO TO 20 301 

BOFEP * B2(P * SORT (E**2 - 1.0) ) 302 

RETURN 303 

20 NOE * NEPS(K) 304 

NOP = NEP ( K ) 305 

I F ( E .GT. 1.0001) GO TO 25 306 

JEPS = 1 307 

23 IQUAL = 2 308 

GO TO 37 309 

C FIND EPSILON BOUNDS FOR INPUT E 310 

25 DO 30 J = 1 » NOE 311 

JEPS = J 312 

I F ( E - EPSILON(J.K) J35.23.30 313 

30 CONTINUE 314 

35 I F ( (K.EQ.2) .AND. (NOE .EQ.JEPS) ) JEPS = JEPS “ 1 315 

C FIND P BOUNDARIES FOR INPUT P 316 

37 DO 40 I = 1 .NOP 317 

IP * I 318 

IF(PTEMP - PF ( I * K ) )50. 45*40 319 

40 CONTINUE 320 

45 GO TO (48.46) .IQUAL 321 

46 BTEMP = BF( IP. JEPS. K) 322 

I F ( K .NE. 2) GO TO 200 323 

BOFEP = BTEMP 324 

RETURN 325 

48 ET ( 1 ) = BF( IP.JEPS-l.K) 326 

ET ( 2 ) = BF( IP.JEPS.K) 327 

IF (K .EQ. 2) ET ( 3 > * BF ( IP * JEPS4-1 *K ) 328 

GO TO (90*72*100) *K 329 

50 I F ( NOP .FQ. IP) IP * IP - 1 330 

GO T 0 160*55) .IQUAL 331 

55 BTEMP = TERP(PF( IP-1. K) ,BF( IP-1.JEPS.K). PTEMP) 332 

IF( K .NE. 2> GO TO 200 333 

BOFEP = BTEMP 334 

RETURN 335 


“k 1 • 



60 60 TO (75*65*75) *K 

C QUADRATIC INTERPOLATION IN P DIRECTION FOR K « 2 
65 DO 70 L * 1*3 

ET(U « TERP<PF(IP-1.IC)*BF<IP-1*JEPS+L-2*K).PTEMP) 

70 CONTINUE 

C QUADRATIC INTERPOLATION IN EPSILON DIRECTION FOR K » 2 

72 BOFEP - TERP(EPSILON(JEPS“l*K) *ET*E) 

RETURN 

C QUADRATIC INTERPOLATION IN PF DIRECTION FOR K » 1 OR 3 

75 DO 80 L ■ 1*2 

ET(L) * TERP(PF(IP-1»K)*BF(IP-1*JEPS+L-2*K)»PTEMP) 

80 CONTINUE _ „ 

C LINEAR INTERPOLATION IN EPSILON DIRECTION ( L0G( 1/ ( E-l) ) FOR K * 1 * 

C LOG(E) FOR K * 3) 

I F ( K .EQ. 3) GO TO 100 

90 TEMP -ALOG (1.0/(EPSILON(JEPS-1*K)-1,0) ) 

BTEMP * E T C 1 ) + ( (ALOG(1.0/(E“1.0) )- TEMP )/( ALOG( 1*0/ 

# (EPSlLONt JEPS.KJ-1.0) )- TEMP)) * (ET<2) - ET<1) ) 

GO TO 200 

100 TEMP * ALOG(EPSILON(JEPS-l.K) ) 

BTEMP * ET ( 1 ) + ((ALOG(E) - TEMP) / ( ALOG ( EPS I LON ( JePS * K ) ) - TEMP 
*) ) * C ET < 2 ) - ET(1) ) 

200 CALL KlK3(E*P»BTEMPl*PTEMPtBTEMP*2*K) 

BOFEP ■ BTEMP 1 

RETURN 

END 


336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 


SUBROUTINE CONV (EPS*P) 

C 

C SETS UP CONVERSION FACTORS FOR TABLES 1 AND 2 

C COMMON /COLD/ P I *EPSQ*SQUARE .ROOT *ACOSINE *PMC *TERMT 

DATA ( P I * 3.141592654) 

EPSQ = EPS**2 
SQUARE = EPSQ - 1.0 
ROOT « SORT (SQUARE) 

ACOSINE = ACOS (1.0/EPS) 

PMC * PI - ACOSINE 

TERMT * (PMC * ( (2.0+EPSQ)/RO0T) + 3.0 > 

RETURN 

END 


362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 
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SUBROUTINE KIK3 <EPS*P*B*ZXI*YBXI*ISWITCH*K) 376 

377 

I SWITCH ■ 1 CONVERT P TO XI OR Z 378 

* 2 CONVERT BXI OR Y TO B 379 

380 

COMMON /COLD/ Pl*EPSQ*SQUARE*ROOT*ACOSINE*PMc*TERMT 381 

IF ( K .EQ. 3) GO TO 25 382 

5 COEF * PI * SQUARE**1.5 / <2.0* EPSQ) 383 

FNEW « EPSQ/ SQUARE 384 

GO TO < 10*20) *ISWITCH 385 

10 ZXI * P * ROOT / FNEW 386 

RETURN 387 

20 B *YBX I * TERMT / (COEF * 5||E*frfc, 308 

RETURN 389 

25 GO TO (30*^0)*fl6triTCH 390 

30 TEMP = ROOT * (Pl/2.0) * SQUARE»*2.5 /<EPSQ**2* TERMT) 391 

ZXI ■ P * TEMP 392 

RETURN 393 

40 TEMP * (4.0*EPSQ*TERMT ) / (PI * SQUARE**1.5> 394 

3 * YBXI * TEMP * B2(P * ROOT) 395 

RETURN 396 

END 397 



FUNCTION BASYM (E.P.K) 398 

399 

CALCULATES ASYMPTOTIC VALUE OF B(P*E) FOR P GT P OF TABLES 400 

SEE PAGE 11* (CM II) 401 

402 

COMMON /COLD/ P I .EPSQ .SQUARE »ROOT *ACOS INE *PMC * TERMT 403 

IFOC.EQ. 2) CALL CONV (E»P> 404 

TB * PMC/ <8.0*ROOT) 405 

TC = 9.0 *EPSQ**3 + 138,0 *EPSQ**2 + 152,0*EPSQ +16.0 406 

TD = ( 1.0/24. 0)*(229.0*EPSQ**2+592.0*EPSQ+124.0) 407 

TA = EPSQ/(2.0*P**3 * SQUARE**5) 408 

TERM1 = (EPSQ/(2.0*P * SQUARE**2>) * TERMT 409 

TERM2 = TA * (TB * TC + TD) 410 

I F ( TERM2 .GT. ( 1.0/3. 0 > *TERM1 ) TERM2 * 0.0 411 

BASYM * TERM1 + TERM2 412 

RETURN 413 

END 414 




r 
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FUNCTION TERP ( X * Y*X IN ) *15 

C QUADRATIC INTERPOLATION *16 

DIMENSION X<3)*Y(3) *17 

SUM * 0.0 *18 

DO 100 J * 1*3 *19 

PROD * 1.0 *20 

DO 50 I *1*3 *21 

IF ( I .EQ. J> GO TO 50 *22 

PROD * PROD * UXIN-X(I) 1 / (X(J) -X(D ) ) *23 

50 CONTINUE *24 

SUM » SUM + PROD * Y ( J ) *25 

100 CONTINUE *26 

TLRP = SUM *27 

RETURN *28 

END *29 


FUNCTION B2<Z) *30 

C *31 

C SEE PAGE 11* (C) (II) 432 

C *33 

DIMENSION ZZ(20) *BB(20) *NJ5) *C(5) *34 

DATA ( ZZ = .0* .1* .2* *4* .6* *8* 1.* 1*2* 1*4* 1*6* 1.8» 2.* 3.* 435 

* A.*5.*6.*7.*8.*9.*l0.) * *36 

* (BB * 0.* 5.8501E-2* 1.6039E-1. 3.5887E-1. 4.9819E-1* *37 

* 5.7597E-1* 6.0588E-1 * 6.0315E-1* 5.8021E-1* *38 

* 5.4590E-1* 5.0663E-1* 4.6692E-1* 3.0540E-1* *39 

* 2* 1580E-1 * 1.6640E-1* 1.3575E-1* 1.1509E-1* **0 

* 1.0005E-1* 8.8559E-2* 7.9470E»2> **1 

DATA (N=l*3*5»7»9) * ( C * .7854 . .8836 ,4. 142 *42. 28,750 .* ) **2 

IF ( Z .GT.10.0 ) GO TO 200 443 

100 DO 125 I = 1*20 **4 

IF ( Z - ZZ(i> ) 110*105*125 445 

105 B 2 * BB( I) 446 

RETURN 447 

110 B2 * BB(I-l) + ( (Z-ZZ( 1-1) )/(ZZ( I )-ZZ( 1-1) ) ) * < BB Cl) - BB<I“1) ) 448 

RETURN 449 

125 CONTINUE 450 

200 B2 = 0.0 451 

DO 225 I * 1*5 452 

B2 * B2 + C ( I ) / Z**N( I ) 453 

225 CONTINUE 454 

RETURN 455 

END 456 




